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Abstract. In the paper, we study entire solutions of the difference equation t/> (z + h) = M (z) ip (z), 
"^j". z £ C, if) (z) £ C 2 . In this equation, h is a fixed positive parameter, and M : C m SL (2, C) is a 

' given matrix function. We assume that M(z) is a 27r-periodic trigonometric polynomial. The main 

aim is to construct the minimal entire solutions, e.i. the solutions with the minimal possible growth 
' simultaneously as for z — > —ioo so for z — > +ioo. 

We show that the monodromy matrices corresponding to the bases made of the minimal solutions are 
trigonometric polynomials of the same order as the matrix M . This property relates the spectral anal- 
ysis of the one dimensional difference Schrodinger equations with the potentials being trigonometric 
C*~) ' polynomials to an analysis of a finite dimcnsinal dynamical system. 

> 

O 1. Introduction: Some definitions and main results 

(N ' 

O ■ The paper is devoted to the study of entire solutions of the equation 

(1.1) ip(z + h) = M(z)4>(z), zeC, f(z)eC 2 . 

In this equation, h is a fixed positive parameter, and M is a given matrix. We assume that the 
matrix M satisfies the following two conditions. First, 

> 



M (z) G SL (2, C), zeC, 
secondly, the matrix M is a 27r-per iodic trigonometric polynomial. 

Since equation ( |1 . 1| ) is invariant with respect to the multiplication of the solutions by the /i-periodic 
functions, one can construct its entire solutions growing as for z — > —ioo so for z — > -Moo as quickly 
as wanted. The main aim of this paper is to construct the minimal entire solutions, e.i. the solutions 
with the minimal possible growth simultaneously as for z — > —ioo so for z — > +ioo. 

The set of solutions of ( |1 . 1|) is invariant with respect to the operator of 27r-translation: f(z) i— > 
f(z + 2Tc), and there are natural objects, monodromy matrices, describing the transformations of the 
bases in the space of solutions of (|1.1|) under such translations, see subsection pL . 10] . The notion of 



X 



the monodromy matrix for difference equations with periodic coefficients is a natural generalization 
of the classical notion of the monodromy matrix for ordinary differential equations with periodic 
coefficients. For the differential equations, the monodromy matirces are constant, i.e. independent 
of the variable of the equation. For the difference equations, the monodromy matrices are periodic 
functions of the variable. We show that the monodromy matrices corresponding to the bases made of 
the minimal solutions have the most simple analytic structure: they are trigonometric polynomials 
of the same order as the matrix M. 

The monodromy matrices appear in the spectral analysis of one dimensinal Schrodinger equations 
with periodic potentials. The spectral analysis of a differential Schrodinger equation with a periodic 
potential reduces to the study of the spectrum of a (constant) monodromy matrix. This leads a simple 
"discret" band structure of the spectra of the periodic differential equations. For difference equations, 
the monodromy matrices being periodic, instead of the spectral analysis of an individual monodromy 
matrix, one arrives to an infinite sequence of the monodromy matrices and has to study properties 



of an infinite sequence of finite difference equations of the form (LI). This reflects the cantorien 



structure of the spectra of difference equations with periodic coefficients. The above property of 
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the monodromy matrices corresponding to the minimal entire solutions is important for the spectral 
analysis of the difference Schrodinger equations with the potentials being trigonometric polynomials. 
It relates their spectral analysis to an analysis of a finite dimensinal dynamical system. 

This paper is a natural continuation of the articles 0, [|J, 0, 0, @] and devoted to Harper 
equation 

/„ ~\ ib (z + h) + ib (z — h) , , x ,_, , , . 

(1.2) — ; r v i + cos 2^0*0 = Eib(z) 

and was inspirated by the papers of B.Helffer and J.Sjdstrand, and of Wilkinson (see, for example, |§, 
||). They have suggested an asymptotic renormalization method to study the spectrum of the Harper 
equation (|1.2j ). Our works are related to an exact renormalization procedure, the monodromization 
method. 

1.1. Notations and agreements. 1. To simplify the notations, we let 

a (z) = Mn(z), b (z) = M 12 {z), c (z) = M 2l {z), d (z) = M 22 (z). 



2. If either b(z) = 0, or c(z) = 0, then the vector equation ( |1 . 1| ) reduces to two scalar equations 
which can be solved explicitly. In this paper, we always assume that 

(1.3) b(z)^0, 
and let 

(1.4) p(z)=b(z)/b(z-h), v(z) = a(z) + p(z)d(z- h). 

These two functions play an essential part in our constructions. Elementary calculations show that 
if ib is a vector solution of ( |1 . 1| ) , then its first component satisfies the equation 

(1.5) ibi{z + h)+p (z) ip x {z -h)=v (z) ^(z). 

3. When describing the vector solutions of ( | . 1|) , we describe only their first components: if the 
first component of a vector solution is known, the second one can be reconstructed by the formula 

(1.6) M*) = ^ Z + V~ a ^^ z \ 

b(z) 

4. Throughout the paper, we use also the following notations. Let / be a 27r-periodic function. If 

/ (z) = f. m e- imz (l + o (1)), /_ m ^0, z - +200, 

we put 

n+(f) = m > /+ = f-m, 

and if 

f{z) = f l e a *(l + o{l)), fi^O, z^-ioo, 

we let 

n-{f) = l, f- = fi- 

1.2. Set of analytic solutions of (|1 . 1| ) . List some elementary properties of the set of analytic 
solutions of equation fll.lj ), see, for example, M. 

Fix a, (3 G R, a < p. Let S be the strip 

S = {z E C : a < Imz < (3} 

and let M (S) be the set of the vector solutions of ( |1 . 1|) analytic in this strip. Denote the ring of all 
/i-periodic functions analytic in S by K (S). 

For any two solutions ibi, ib 2 G M (S) the wronskian 



is an element of K(5 r ). The solutions ?/>i, ip2 are linearly independent over the ring K (S) if 
Wz),^)}/0 )2 g5, 

Let two solutions ipi, ip2 £ M (S 1 ) be linearly independent over the ring IK(5'). Then any solution 
i/)6M (S) can be uniquely represented in the form 

(1.8) V = 0^1 + ^2, a,peK(S). 

Thus, the set M (S) is a two dimensional module over the ring K (S), and the solutions tpi, ip2 form 
a basis of M(S). 

In the case of S = C, we write simply M and K. 

1.3. Bloch solutions. The set M (S) is invariant with respect to the 27r-translations. We call 
/ G M (S) a Bloch solution if 

(1.9) f(z + 2ir)=u(z)f(z), ueK(S). 

The factor -u is called the Bloch multiplier of the solution /. Bloch solutions play an important role 
in the theory of (|1.1|) . 

1.4. Singular points of equation (|1 . 1| ) . The 27r-periodicity of the matrix M makes natural to 
consider -Moo and —zoo as two singular "points" of the equation. For any fixed Y e M, call the 
half-plane C + = C + (Y) = {zgC : Imz > Y} a vicinity of -Moo, and the half-plane C_ = C-(Y) = 
{z E C : Imz < Y} a vicinity of the —zoo. For brevity, we shall write 

M± = M(C±), K± = K(C±). 

The minimal entire solutions of equation ( |1 . 1|) can be characterized as the entire solutions having 
the "simplest" asymptotic behavior as for z — > —ioo so for z — > -Moo. 

1.5. Canonical bases in vicinities of the singular points. We shall define the minimal solutions 
in terms of the "simplest" solutions living and analytic in some vicinities of the singular points ±zoo. 
They appear to be Bloch solutions. One has 

Theorem 1.1a. Let b ^ 0, and let v — > oo as z — > -Moo. There exist two vector Bloch solutions f\ t 2 
of ( |1 . 1|) analytic in a vicinity C + of -Moo with the first components admitting the representations 

(1.10) {h 2 \{z) = e ± ^M (n+ W z + ^ + % ( n+ ^ ~ t + 
as z — ► +ioo. Here, 

(1.11) (j) + = i\nv + -±n + {b). 
These two solutions are linearly independent over the ring K +; 

(1.12) det(A(*),/ 2 (z)) = ^. 

Theorem 1.1b. Let b ^ 0, and let v — > oo as z — > —ioo. There exist two Bloch solutions of ( [L.l|) 
analytic in a vicinity C_ of —ioo with the first components admitting the representations 

(1.13) {gi 2)^) = e ± ^M ^ v)z + ^ ~ l ^ v) ~ f + W, 
as z — > —ioo. Here, 

(1.14) = -ilnu_ - |n_(6). 
These two solutions are linearly independent over the ring K_, 

(1.15) det(^(z),(/ 2 (z)) = -^. 



Remark 1. Note that the numbers <j>± are defined only modulo 2ir. Fixing, for example, the parameter 
in two different ways, we obtain two different pairs fx 2 . 

Remark 2. The asymptotic formulas for fx 2 and g± immediately imply the asymptotics of the corre- 
sponding Bloch multipliers. For /i )2 the Bloch multipliers admit the representations 

(1.16) ai, 2 (z) = al 2 e ±27rn +(^ iz / h (l + o(l)), z - +ioo, 
and, for gi i2 , the Bloch multipliers have the form 

(1.17) /? 1 , 2 (z)=/? 1 ° 2 e ±2 ^-(^/ /l (l + o(l)), z^+ioo, 
here a\ 2 an d 0x2 are nonzero constants, 

"1,2 = ex P 
Pi,2 = exp 



27TZ 

±_ /T ^ + + 7rn +( w )) + i7r ( n +( v ) ~ n +( & )) 



27TZ 

±-^— (0_ + nn~(v)) — i7r(n_(u) — n_(b)) ) . 



We see that the pair / 12 is a basis of M + , and the solutions gi j2 form a basis of ML. We call these 
bases canonical. We call the number + the parameter of the canonical basis /i i2 and the number 0_ 
the parameter of the canonical basis #i i2 . 

1.6. Consistent canonical bases. Let 

S = {±(2<K + h + 2hl + 2Trm), l,m E N}. 
We call two canonical bases /i j2 and #i j2 consistent if their parameters satisfy the condition 

(1.18) 

Since the numbers 4>± are denned modulo 2tt, the consistent bases always exist. 

1.7. Minimal entire solutions. Existence. Let the assumptions of Theorems la - lb be fulfilled. 

Fix two canonical bases fx t 2 and gx t 2 analytic in some vicinities C± of ±200. Any entire solution of 
fll.lj) admits the representations 

(1.19) 4>(z) = A(z)fx(z)+B(z)f 2 (z), zeC + , A, B e K + , 

(1.20) ^ (*) = (7 (s)0i (*)+!> (3)02(3), z6C_, C.flGL. 

Note that the periodic coefficients in these representations are uniquely defined. 

We call a nonzero entire solution of ( |1 . 1|) minimal, if 
( i) its coefficients A, B are bounded in C+; 
( it) the coefficients C and D are bounded in C_; 

(Hi) at least, one of these four coefficients tends to zero as Imz tends to the corresponding 
infinity. 

Note that this definition depends on the choice of the canonical bases. For a given pair of canonical 
bases, the definition singles out at least four minimal solutions. 

One of the main results of the paper is 
Theorem 1.2. Let b(z) ^ 0, and 

(1.21) n + (v) = n-(v) = n > 0. 

Then, for any two consistent canonical bases, there exist all the four corresvondinq minimal entire 



Note that under the assumptions of this theorem, there exist all the canonical bases. We have already 
explained the condition on the coefficient b. Discuss the assumption ( [L.21 ). We shall see that among 
the equations ( |1 . 1| ) , there is a nontrivial one which can be solved in terms of certain contour integrals. 
This is the equation with 

(1.22) a(z) = -2A cos{nz), b(z) = -1, c{z) = 1, d(z) = 0, 

where n G N, and A is a complex number. The proof of the existence theorem is based on this 
observation. In the case of ( |1.22| ), v (z) = 2\cos(nz). This leads to ( |1.21| ). 

Now, we shall discuss the basic properties of the minimal solutions. 

1.8. Asymptotic coefficients of the minimal solutions. Consider the representations (|1.19| ) for 
a minimal solution. The coefficients A, B, C and D can be represented by the converging Fourier 
series 

oo oo 

A (z) = Ae 2mlz / h , B{z)=Y, B ie 2mlz / h , z e C+, 
z=o z=o 



oo oo 



C(z) = J2 Qe- 2mlz / h , D(z)=J2 D ie - 2mlz / h , z e C_. 

1=0 1=0 

One of the coefficients A , B Q , Co and D is zero. 

Assume that the Fourier coefficient Aq equals to 0. In this case, we denote the minimal solution by 
if) a, and call the Fourier coefficients 

Ai, B , Co and D 
the asymptotics coefficients of this minimal solution. 

If Aq 7^ 0, but Bo = 0, we denote the minimal solution by ipB, and call the Fourier coefficients 

A Q , Bi, C and D 

the asymptotics coefficients. 

Continuing in the same way, we define the minimal solutions ipc an d ipDi and their asymptotic 
coefficients. 

In the sequel, we denote the asymptotic coefficients of the minimal solution ipA by A a, Ba, Ca, 
Da, the asymptotic coefficients of the minimal solution ipB by Ab, Bb, Cb, Db and so on. 

The asymptotic coefficients play a crucial part in the analysis of the minimal entire solutions. Let 
us describe some elementary observations. 

1.9. Wronskians of the minimal solutions. Let there exist two canonical bases and in 
some vicinities of ±ioo. Denote the wronskian of /12 by Wf and let w g be the wronskian of g\^- 
Assume also that there exist all the corresponding minimal solutions. In Section 6.1, we prove 

Proposition 1.3. The wronskians of any two of the minimal solutions are constant. One has 

{ipA, ^Pb} = -w f B A A B , 
{ip A , ipc} = —Wf Ba Ac = w g C A D c , {if>A, ^d} = -w f B A A D = -w g D A C D , 
Wb, M = w f A B B c = w g C B D c , {ip B , iPd] = Wf A B B D = -w g D B C D , 

We, *Pd} = -w g D c C D . 



Corollary 1.4. If all the asymptotic coefficients are nonzero, then any two of the minimal solutions 
form a basis o/M. 

In Section 6.2, we obtain also 

Corollary 1.5. (Uniqueness Theorem) If all the asymptotic coefficients are nonzero, then any 
of the minimal solutions is unique up to a constant factor. 

The existence theorem is the central result of this paper. Remind that the matrix M is a trigonometric 
polynomial. In the next publication, we shall study a family of the equations ( |1 . 1| ) parametrized by 
the constant coefficients of M. In particular, we shall prove that the minimal solutions and their 
asymptotic coefficients can be considered as meromorphic functions of these parameters, and that 
these functions are not identically zero. Thus, we shall see that the asymptotic coefficients are non 
zero for any typical matrix M. 

1.10. Monodromy matrices. The notion of the minimal solution is the first of the main notions 
of this paper, and the second one is the notion of a monodromy matrix. We begin by recalling the 
general definition of a monodromy matrix and the description of the monodromization procedure, 
see||. 

Now, we do not need to assume that the matrix M is a trigonometric polynomial. In fact, we have 
to suppose only that it belongs to SL (2, C) and is a 27r-periodic matrix function. 

1. Let \1/ is a matrix solution of (|1 . 1|) . We call this solution fundamental if det \l/ (z) = Const ^ 0. 
Note that ^(z + 2tt) is a solution of ( |1 . 1|) together with ^{z). We define the monodromy matrix 
corresponding to a given fundamental solution \1/ by the relation 

^(z + 2tt) = ty(z)M{(z), 

where * is the transposition. The function Mi has the properties: 

(1.23) detMi(z) = l, M 1 (z + h)=Mi(z). 

Let ip^ and ip^ be two vector solutions of ( |Llp . Compose of them the matrix ?/> < - 2 ' ) ). If 

its determinant is a nonzero constant, then this matrix is a fundamental solution. We define the 
monodromy matrix for such ipi and ip2 as the monodromy matrix corresponding to the fundamental 
solution ip^). 

The notion of a monodromy matrix is well known in the theory of the ordinary differential equations 
with periodic coefficients. For a differential equation 

=M(z)#(z) 

with a 27r-periodic matrix M, the monodromy matrix also defined by ^ (z + 2-k) = (z) M*, but, it 
is independent of z. 

2. Recall the description of the monodromization procedure. Consider the sequence of the numbers 
hj, j — 0,1, ... , defined by the relations 

ho = 2tt, hi = h, 

hj-i = pjhj + hj+i, < hj + i < hj, pj G N. 
The pj are the denominators of the continued fraction 

, „ 1 



Let \l/o be a fundamental solution of equation (Q|), and let Mi be the corresponding monodromy 



matrix. Bring into the consideration the equation 

(1.24) $ 1 (z + h 2 )=M l (z)V 1 (z). 

In view of (|1.23|) , this equation is of the same type as 



The passage from (|1.1|) to ( |1.24|) is the first step of the monodromization procedure. If there exists 



a fundamental solution of (|1.24|) , the monodromization procedure can be continued. In result, one 
obtains the suite of the monodromy matrices satisfying the relations: 

det Mj{z) = 1, Mj(z + hj) = Mj(z). 

If the number hi /ho is rational, the procedure is finite, In general case, the monodromization procedure 
is infinite. 

The spectral analysis of a differential equation with periodic coefficients reduces to the study the 
spectrum of a (constant) monodromy matrix. Trying develop similar ideas for the difference equa- 
tions, arrives to the infinite sequence of finite difference equations ( |1 . 23|) . In @]-[[|, we have seen 



that spectral properties of (|1 . 1|) with respect to the parameters of M are tensely related to certain 



properties of the sequence Mi, M2, M3, etc. 

1.11. Family of matrices invariant with respect to the monodromization. Denote by r m i , 
m, 1 G N, the set of the trigonometric polynomials of the form 

/ (z) = e~ imz f- m + e~ l[m ~ ^ z f- m+1 + ■■■ + " ^/j-i + e ilz f h 
e.i. such that < I, n + (f) < m. 

Fix nGN. By fl (n) or, briefly, by fl, denote the set of matrix functions M (z) G SL (2, C) such that 

o G T n ri , b G 7Vi,ti— l> 

C G Tn—\,ni d G Tfi— X,n— lj 

and 

n+(a) = n_(a) = n. 



Note that the conditions (|1.3|) and ( 1.21|) of the existence theorem are fulfilled for any matrix from 
the family Q. 

One has 

Theorem 1.6. Let M G fl(n). Consider two canonical bases and g\.2- Assume that the 
corresponding minimal entire solutions ipo and tps exist and that their wronskian is nonzero. Then 
the corresponding monodromy matrix as a function of the variable z\ = 2nz/h also belongs to fl (n). 

This theorem can be used as the base for a program for the spectral investigation of the difference 
Schrodiner equations with the trigonometric polynomial potentials: the equation 

ip(z + h)+ip(z-h) 

2 +p{z)ip {z) = Eip (z) 

(2E — 2v ( z) —l\ 
Y q J • I n this case, M G £1 (n), and the 

monodromization procedure leads to a sequence of equations ( |1 . 1| ) with the matrices from the family 
fl (n). And thus, it is equivalent to a finite dimensional dynamical system. In this paper, we do not 
consider any spectral problems, and note only that the monodromization procedure is close to the 
asvmptotic renormalization approach suggested for ftL.21) by Helffer and Siostrand Hi and Wilkinson 



We finish this discussion by listing the asymptotic formulae for the coefficients of the monodromy 
matrix from Theorem 1.6: 



(1.25) 
(1.26) 
(1.27) 

(1.28) 

(1.29) 

(1.30) 

(1.31) 

(1.32) 



M n {z) 
M n (z) 
M 12 {z) 

M 12 {z) 

M 21 (z) 

M 2 i(z) 

M 22 (z) 

M 22 (z) 



a° 2 e- 2mnz / h (l + o{l)), 
(3° 1 e +2mnz / h (l + o(l)), 
_ a o e -2mnz/h (^ + (1) 

_0O e +2m (n - 1) z/h (d± + Q ^\ 



z 
z 
z 



a o e -2m (n 



1) Z/h ( e AnH/h^ + o(r 



o fi — 2-ni (n — 1) z/h ( ^AiPilh ApB B i Q 



•a, e 



-27ri (n 



1);:/ft (§§§f + °(i)), 



+ioo, 
— zoo, 
+ioo, 

—zoo, 

+ioo, 

—zoo, 

+ioo, 

— ioo. 



Here a? 2 and /3° 2 are the constants from the asymptotic representations (|1.16|) - ( [L.171 ) for the Bloch 
multipliers of the solutions fi j2 and gi t2 . Note also that, since the wronskian of ipo and ipB is nonzero, 
then, in view of Proposition 1.3, Ab, Db, Bd, Cd ^ 0. 

1.12. Typical properties of the minimal solutions. Above, we have described the main results 
of this paper. We shall continue it in the next publication. There, we shall study some typical 
properties of the minimal solutions and, in particular, of their asymptotic coefficients. In fact, the 
second part will be devoted to the study of the dependence of the minimal solutions and their 
asymptotic coefficients on the parameters of the matrix M. Let us formulate here the central results 
of this second part. 



Together with the set fl (n), we consider also its subsets fl m i, 
the conditions 

n + (b) = m, n_{b) = I. 



-n < —m < I < n — 1, singled out by 



Clearly, 



n 



•ml 



ml i 



and 



Qml n 



(m,l)^(m',l'). 



It appears that the elements of Q can be uniquely parametrized by the constant coefficients of the 
trigonometric polynomials a (z) and b (z). Denote by to the set of constant coefficients of the matrix 
M. One can consider u as a point in C 8 ™. With this, the set f2 m /, — n < — m < I < n — 1, appears 
to be a connected analytic submanifold of C 8n of the dimension 2n + m + / + 2, and the set Q (n) 
itself is a connected analytic submanifold of C 8n of the dimension 4n + 1. 

Let D be a simply connected domain of Q m i. Fix some continuous on T> branches of the functions 

<p_ = — I h/2 — i lna +n , <j) + = —mh/2 + i In a_ n . 

Assume that, for any point u G V, these functions satisfy the condition ( |1.18| ), e.i. that the canonical 
bases with these parameters are consistent for all to G T>. One has 



Theorem 1.7. For the above set T> and the above canonical bases, the minimal entire solutions can 
be normalized so that they would be meromorphic inu eD together with their asymptotic coefficients. 
In this case, the asymptotic coefficients are not identically zero. 



Corollary 1.8. In the case of Theorem 1.7, the monodromy matrix corresponding to the minimal 
solutions ipD o,nd ipB is meromorphic in uj GD. Moreover, this matrix typically belongs to fi njn _i. 

Here, we have used the usual terminology of the analytic set theory: one says that a property of a 
function defined on an analytic set is typical if it takes place outside some analytic subset of a smaller 
dimension. 

1.13. The plan of the paper. In section 2, we prove Theorem 1.1. Also, we study some additional 
properties of the Bloch solutions. 

Section 3 is devoted to the scalar equation 

(1.33) fj, (z + h) + /j, (z — h) + 2A cos z n (z) — 0, zGC, 

where h > and A G C are two fixed parameters. We construct and study minimal entire solutions 
of this equation. The solutions are explicitly described in terms of certain contour integrals. Note 
that the case of A = — 1 was already treated in 0. 

Section 4 is devoted to the analysis of the equation 

(1.34) f(z + h) + f(z-h) +2\coszf(z) =w(z)f(z), z G C, 
where w is a meromorphic function satisfying the estimates 

\W(Z)\ < Const e -(l-*)|Im Z | e> n 



for sufficiently big |Im^|. In terms of the minimal solutions of ( |1.33|) , we invert the operator in the 



left-hand side of ( |1.34| ), and obtain an integral Fredholm equation on a contour in the complex plane. 



This allows to construct and investigate some meromorphic solutions of ( |1.34| ) 



In Section 5, we reduce the matrix equation ( |1 . 1| ) to the form ( |1.34j ), and then prove Theorem 
1.2. 

In the first part of Section 6, we study basic properties of the minimal entire solutions assuming 
that their asymptotic coefficients are nonzero. In particular, we prove their linear independence 
over the field of /i-periodic function, their uniqueness (up to independent of z factors), and justify 
Theorem 1.6. Recall that the monodromy matrices corresponding to the minimal solutions ipo and 
ipB are trigonometric polynomials of 2nz/h. Formulae ( |1.25| ) - ( |1.32j ) allow us to calculate only some 



of the constant coefficients of these polynomials. In the end of the section, we describe a way to 
calculate all the other coefficients. 

Section 7 is devoted to some remarks for the case where M is a trigonometric polynomial of the 
first order. 
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2. Bloch solutions in a vicinity of the infinity 

In this section, we construct Bloch solutions / 1)2 of equation ( |1 . 1| ) described in Theorem 1.1a. 
Subsections 2.1 - 2.5 are devoted to the proof of this theorem. 

In subsection 2.6, we discuss the set of all the solutions analytic in C + (C_) and having there 
the same asymptotic representations as /i )2 (#1,2) • 

2.1. The plan of the proof. Consider equation (|1.5|) for the first component ipi of a vector solution 
■ip of equation ( |1.1| ). Setting 

(2.1) $(z) = M* + h)/M*)> 

we see that 



We call (|2.2|) a difference Ricatti equation. 

To construct a Bloch solution of ( |1 . 1| ) analytic in a vicinity of +ioo we use the ideas of ||. 
The first step is to prove the existence of an analytic 27r-periodic solution $ of the difference Ricatti 
equation in a vicinity of -Moo. Then we introduce the function 

(2.3) <j>(z)=\n<S>(z) 
and consider the equation 

(2.4) \(z + h)-\(z) = <j)(z). 

Having solved this equation we reconstruct the first component of a vector solution of ( |1.1[ ) by the 
formula 



(2.5) ^(z) 



after that the second one can be recovered by ( |1.6| ). Finally, we check that the 27r-periodicity of $ 
implies that the constructed solution is a Bloch solution. 

2.2. Analytic solutions of the Ricatti equation. Proposition 2.1. In the case of Theorem 
1.1a, there exist two Ik -"periodic solutions $1,2 of (|2.2|) analytic in a vicinity of +ioo and having the 
asymptotics: 

(2.6) $ 2 = v(z) + o(l), z^+ioo, 
and 

(2.7) ^(z) = ^f±^-(l + o(l)), z^+ioo. 

v [z + h) 

Proof. Consider the sequence of the functions 

&°\z)=v(z), *<»{z)=v(z)- , P(Z) , V 

v(z — h) 



The limit of this sequence (if exists) is a solution of ( |2.2|) represented by the continuous fraction 
(2.8) * a (z)=v(z) f,{:) 



v (z — h) P( z ~ h ) 



This continuous fraction was investigated in |J. 

Recall that the coefficients of the matrix M are trigonometric polynomials, and that 6^0. Both the 
functions v and p are analytic in a vicinity of +ioo. Moreover, for z +ioo, the function p tends 
to a finite limit, and v (z) — > oo. Let 

v-(Y)= inf \v(z)\, p+(Y)= sup \p(z)\, (i = P+ 

lmz>Y lmz>Y l V -/ 2 ) 

As in 0, one can easily check that if p, < 1, then 

- $W( Z )| < // |$W(z) - ( z ) |, Imz>7. 

This estimate implies that all the functions n = 1, 2, . . . , are analytic near the "point" +ioo, 

that the continuous fraction in (|2.8| ) converges uniformly in z in the half-plane C + (Y) for sufficiently 
big Y, and that <£>o satisfies the estimate (12.61). Since anv of the functions $ ( - n ^ is 27r-periodic, the 



The solution $1 can be represented by the continuous fraction 

p(z + h) 



v(z + h) V(z + W) 



It can be investigated in the same way as Q2.8| ). □ 

2.3. Logarithms of the solutions of the Ricatti equation. Let 

(2.9) 0i, 2 (z) = ln$i i2 (z), 

where <£>i i2 are the solutions of the Ricatti equation constructed in Proposition 2.1. Asymptotic 
representations fl2.6Q - (|2.7|) immediately imply 



Lemma 2.2. In a vicinity of+ioo, the functions <fiifl{z) are analytic and can be represented in the 
form 

(2.10) 01,2(2;) = ±in + (v) z + a 1)2 + 51,2(2), 
where 

(2.11) a 2 = In v + , ai = — In v + — in + (b) h + in + {y) h, 
and g\ t 2 are 2tt -periodic analytic functions decaying as z — > +ioo. 

2.4. Homological equation. Here, we collect some facts on the homological equation 

(2.12) (p(z + h)-<p(z)=g(z), z e C, 
where g is a 27r-periodic function. 

1. First, prove 

Lemma 2.3. Let g be a 2ir-periodic function analytic in a vicinity of +ioo decaying as z — > +ioc, 
then equation ( p. 12 ) has a solution analytic in a vicinity of+ioo and decaying as Imz — > +00 
uniformly in Kez if |Rez| is bounded by a constant. 

Proof. The solution (p can be constructed by explicit formulas. Fix a point z Q in the domain of 
analyticity of the function g. Let 

1 rz +2-K 

G(k) = — e- lkz g(z)dz, 

and 

(2.13) v(z)= l.jG {k)e ^^± dk , 
Here 7 is a contour in the complex plane of k. Describe it. 

The integrand in ( |2.13| ) has poles only at the real line. One of the poles is at k = 0. One is situated at 
the point k Q = min{l, 2ir/h}. All the other are located outside the interval [0, k }. Fix two constants 
Ci and c 2 so that c\ > 0, < c 2 < k Q . The contour 7 is a contour coming from +00 along the line 
Imz = Ci , going from the upper half-plane of the complex plane along the line Re z = c 2 , and coming 
back to +00 along the line Imz = —C\. 

If Im z > Im z Q , the integral for ip converges absolutely. Moreover, 

\<p(z)\< Const e -c 2 Imz_ 



One checks that <p satisfies equation (|2.12|) by means of the residue theorem: 

oo 

<p(z + h)-(p(z) = — G{k) e lkz ctg (irk) dk = J^ e "" G (k) = g (*)■ 

2l A K=1 

□ 

2. Note that if <p is a solution of ( |2.12| ), then the function g(z) = p> (z + 2tt) — <p (z) is /i-periodic: 

g (z + h) — g (z) = <p(z + 2ir + h) — <p(z + h) — <p(z + 2tt) + p> (z) = 
= g(z + 2n)-g(z) = 
= 0. 

2.5. Bloch solutions in a vicinity of -Moo. Now, we have all the ingredients to construct two 
Bloch solutions of (|1.1| ) analytic in a vicinity of the point -Moo. 



2.5.1. First solution. Let $2 be the solution of the Ricatti equation fl2.2|) described by Proposition 
2.1 and let <j) 2 be the logarithm of $2- Consider the equation 

(2.14) \{z + h)-\{z) = <p 2 {z). 

Remind that <p 2 admits the representation ( |2.10| ). Therefore, any solution of (|2.14j) can be represented 
in the form 

2 

z z z 

(2.15) \{z) = -in + (v) — + in+(v)- + \nv + - + <p{z), 

where p> is a solution of the homological equation with g 2 (z) in the right-hand side. We construct 
a solution *p of this equation as in Lemma 2.3, and then, we construct a solution ip (z) of ( |1 . 1| ) by 
formulae Q2.5| ) and ( |l.(j| ). The result can be written in the form 

(2.16) ^(z) =e x ^ ( $2(z )- a ( 2 ) 

V b(z) 

Note that since b (z) is a trigonometric polynomial, the function is analytic in a vicinity of -Moo, 
and, therefore, ip is also analytic in some vicinity C + of -Moo. 

Check that ip is a Bloch solution. Since the vector in the right hand side of ([2.1 6|) is 27r-periodic, 

%j)(z + 2vr) = 7(«)V(z), zeC+, 

where 

7(z) =e X(z + 2n)-X(z)_ 

Thus, it suffices to show that 

(2.17) 1 {z + h)= 1 {z). 



By (EH), 



Const e -n + {v)^ e ^z + 2ir)-^z) 



7 \ z ) 

But, in view of section 2.3.2, the function <p(z+2ir) — ip(z) is /i-periodic. This implies the /i-periodicity 
of 7. 

Denote the constructed solution by f 2 . We have come to 

Lemma 2.4. There exists a Bloch solution f 2 of ( |1 . 1| ) analytic in a vicinity of +ioo with the first 
component admitting the representation: 



2.5.2. Second solution. Starting with $1 and repeating the above arguments we obtain 



Lemma 2.5. There exists a Block solution fx of ( |1.1|) analytic in a vicinity of +ioo with the first 
component admitting the representation: 

2.5.3. Linear independence. Calculate the determinant of the matrix F (z) = (fx(z),f 2 (z)). Take 
advantage of formula ( |1.8| ). This gives 

detF(z)= (/l)l( ^ ( ( a)l(z) (^-^(z)). 
&(*) 

Substituting in this formula the asymptotic representations for / 1)2 and $1,2, we come to the repre- 
sentation 

(2.18) detF(z) =v+/b+(l + g(z)), g (z) -> 0, 2 -> +zoo. 

Here g is an /i-periodic function, analytic in a vicinity of +200. Note that formula ( [2.18 ) implies that 



fi and f 2 are linearly independent over the field of /i-periodic functions in a vicinity of +ioo. 

Now, we redefine the Bloch solutions fi and f 2 . 

fi:=qfi/(l + g), / 2 :=/ 2 /g, ? = exp(^ 2 + / (2/m+(t;))). 

The new fx and / 2 are also Bloch solutions of equation ( |1 . 1| ) , but their wronskian is constant. They 
are the solutions described in Theorem 1.1a. The proof of Theorem 1.1b is absolutely similar. □ 

2.6. Uniqueness of the Bloch solutions. 1. First, we answer the question:" Are the Bloch 
solutions described in Theorem 1.1 uniquely defined by their asymptotics for z — > ±ioo ?" 

Lemma 2.6. Let f be a Bloch solution of ( |1 . 1|) analytic in some vicinity C + of +ioo and, as the 
solution fx, having the asymptotics ( l.lOj) for z — > +ioo. Then 

f(z)=c(z)f 1 (z), 

where c(z) is an h-periodic function analytic in a vicinity of +ioo, and having the asymptotics 

c (z) — > 1, z — > +ioo. 

Proof. Since the solutions /i i2 are linearly independent in some vicinity of +ioo, / can be represented 
by their linear combination with some ^,-periodic analytic coefficients 

f(z) = a(z)fx(z)+(3(z)f 2 (z). 

The coefficients a and (3 can be expressed explicitly in terms of the solutions: 

_ jet (/, f 2 ) det(fx,f) 
01 det (fx, f 2 y P det (fx, f 2 y 

Since both fx and / are described by (|1.10Q as z — ► +ioo, 

a (z) — > 1, z — > +ioo. 

Now, show that (3 = 0. Since / and fx are Bloch solutions 

N-l 

Hpt. ff, f\(y. 4- 9.7r/\n = Hpt ff, fW*\ TT 7/.,/'* 4- 9-rr/W/ (z 4- 9tt/^ V /V <= M 



where U\ and u are the Bloch coefficients of f\ and /. The asymptotics ( |f .f 0| ) (see also Remark 2 to 
Theoremsl.l) imply that, for sufficiently big Imz, 

\ui(z)\, \u(z)\ < Ce- 27rn+{v)lmz/h . 



This formula is uniform in Re z since ( |f.i 0| ) is uniform in Kez if \Kez\ is bounded by a constant, 



and since u\ and u are /i-periodic. In result, for sufficiently big Imz, 

| det (fx, f) (z + 2ttA^)| < Ce-^ Nn+{v)lmz/h \ det (fx, f) (z)\, ViVeN, 

uniformly in Re z. Since the determinant is /i-periodic, this is possible only if it equals to 0. Thus, 
f(z) = a (z) fi(z) which proves the lemma. □ 

We have investigated the uniqueness of the solution fx. For the solutions f2, g\ and g 2 one can prove 
the similar statements. 



2. In constructing the Bloch solutions, we have fixed the parameters <j)± defined by formulae ( |I . 1 1 ) 



and (|I.f 4j) only modulo 2n. The choice of the parameter <p + influences the asymptotics of 

fi.2 (91,2) as z —>■ +ioo (z — * — ioo). Fixing this parameter in two different ways, we obtain two 
different bases /i >2 (91,2)- Let us study relations between these bases. We shall indicate explicitly 
the dependence of /i )2 on <p + . On has 

Lemma 2.7. 

(2.19) / 1)2 (2,0 + + 2tt) = c li2 e ±27r ^/ /i / 1 , 2 (z,0 + ), 

where ci )2 = ci j2 (z, 0+) are some analytic h-periodic functions tending to nonzero constants as z — > 
+ioo. 

The proof of this lemma is similar to the proof of the previous one. Again, one can prove the same 
statement for the Bloch solutions gi. 2 . 

The most interesting case is one where we change <p + by 2irn + (v). In this case, one comes to 
Corollary 2.8. 

(2.20) fx, 2 (z, <P + + 2-Kn+(v)) = C^f^z + 2tt , <j>+), 

where Ci >2 = Cx,%(z, <f)+) are some analytic h-periodic functions with the asymptotics 

C 1>2 (z, + ) = e -^K(^) - n + (b)) + Q z ^ +lOQ _ 



3. Minimal solution of the auxiliary equation 
The next three sections are devoted to constructing of the minimal solutions of equation (|1 . 1|) . 



Let us agree on the terminology and the notations we shall use there. 
For z G C, we call the set 

{CGC : Im2 = ImC, -6 < Re (z - () < 5} 

the horizontal 5-vicinity of z. We call the horizontal 5-vicinity of a curve 7 C C the union of the 
horizontal 5-vicinities of the all its points. We call the horizontal distance between a curve 7 and a 
point z the maximal 5 for which the horizontal 5-vicinity of 7 does not contain z. 

We call a curve 7 C C vertical, if it intersects all the lines Imz = Const only at nonzero angles. 
If all these angles are strictly bigger than some fixed positive constant, we call the curve strictly 



Furthermore, in the sequel, we use the letter C as a symbol denoting constants independent of 

z. 

In Section 3, we construct and investigate minimal entire solutions of the model equation 

, . ( m x (z + h)\ _ ( -2etcmz -l\fm 1 (z)\ 

{m 2 (z + h))-{ 1 J\m 2 (z)J> ZE ^ 

where £ is a complex number. Note that the first component mi of its vector solutions satisfies the 
equation 

(3.2) m (z + h) + m (z — h) + 2e^ cos z m(z) = 0, 

and the second component is related to mi by the formula 771-2(2) = 7711(2 — h). The main results 
are formulated in Theorem 3.2. When describing the results, we use the special function described 
in Section 3.2. 

In [0], we have considered the case where £ = 0. Here, we discuss in details only the new 
elements of the proof. 



3.1. Reduction to a first order difference equation. We shall construct a solution m(z) of ( jOj 
in the form: 

(3.3) m (z) = e - iz2 / 2h J e ~ iz P/ h ~ + ™P/ 2h v (p) dp, 

where T is a contour in the complex plane of p. Substituting (|3.3j ) into (j3.2|) , we formally get 

-iz I -izp/h - tp 2 /Ah + mp/2h ^-ip-ih/2 + ^ v (p) dp+ 

-izp/h - ip 2 /Ah + mp/2h (jp-ih/2 + e n v ^ dp = 0; 

or, after having changed the variables, 

e -izp/h - i(p - hf/Ah + m(p - h)/2h r e -i P +ih/2 + e £\ v ^ _ ^ dp+ 

e -izp/h - i(p + hf/Ah + m(p + h)/2h r £ ip+ih/2 + e ^ v ( p + h)dp = 0, 



r 



e +lz I e 



r 



r+h 



r-h 



where T ± h are the contours obtained of T by ±/Vtranslations. Assuming that v is analytic in a 
sufficiently large vicinity of T, and that there are no problems of convergence of the integrals, we 
deform the integration contours back to T and obtain 



e 

r 



-izp/h-ip*/4h+nip/2h ^ e -ip/2+ih/2 + e €+i P /2) v(j) ( e ip/2+ih/2 + ^-ip/2^ + dp = Q 

So, if v satisfies the first order difference equation 

(3.4) v (p + h) — q (p) v (p — h), 

with 

e ih/2 + e ip + i 



q(p) 



; £ 4. e ih/2 + ip' 



then m is a solution of ( |3.2|) . Now, our aim is to construct a solution of Q3.4|) . To get a detailed 
information about this solution we shall need a special function playing an important role in the 



3.2. (T-function. Here, we construct a meromorphic solution of the equation 
(3.5) a(z + h) = (l + e~ iz )a(z-h), z E C. 

A similar function was introduced and systematically used in the diffraction theory ||12j| . Later, it 



was also introduced and studied in a different context by other authors, see, for example, |T3 and 



14]| . So, we discuss this function omitting most of the proofs. 

1. Let C„- be the complex plane cut along the real line from — oo to —tt, and from tt to +oo. Fix 
on C n a branch of the function l (z) = In (1 + e~~ %z ) by the condition 

(3.6) lo( z ) ~ * 0) z ~ *• —ioo. 

Define z 

L Q {z) = / l (z')dz r , z e C^, 

J —ioo 

where the integration contour belongs to C^. Let 
(3-7) fl °(*) = iW TtV^- 



Here, 7 C C n is a strait line Re z = Const passing between the points z ± h. The function # is 
analytic in the strip So = {z e C : — 7r — /i<Rez<7r + /i}. 

By means of the residue theorem, one can easily check that 

9 (z + h) - 6 (z - h) = l Q (z), z±heS . 

Therefore, the function 

a(z) = e 6 ^ 

is a solution of ( |3.5| ) in the strip Sq. By means of ( |3.5| ), one can continue it meromorphically in 
the whole complex plane. Let now a denote this meromorphic solution. Investigate its analytic 
properties. 

2. Discuss the set of poles and zeros of the function a . By construction, the function a is analytic 
and has no zeros in So- In result of the continuation, each zero z = it + 2nl, I = 0, 1, 2, . . . , of the 
function 1 + e~' tz produces a chain of zeros of o at the points 

tt + h + 2ttZ + 2hk, fc = 0,1,2,..., 

and each its zero z = —71 — 2nl, I = 0, 1, 2, . . . , produces a chain of poles of a at the points 

-7r - h - 2-kI - 2hk, k = 0, 1, 2, .... 

3. The pole z = —n — h is simple. One can calculate explicitly the residue of a. Omitting long 
elementary calculations, we write down the result 

/fl in 2 in_ %h_ 
-e 12ft 4 12. 
7T 

One can also get the following explicit formula 

(3.9) a(-n) = ±e^ 2 / ™ + ih/U 

v2 

4. Let us describe the asymptotics of o for Imz — > ±z"oo. Clearly, 

Lq(z) = O (e — 1^1), Imz — > —ioo. 

On the other hand, 



This implies the asymptotics 

(3.10) a{z) = l + o{e~^ lmzl ), Imz^-zoo, 
and 

(3.11) a(z) = e - % + z Tk +i n(l + (e^ |Imz| )), Im^+ioo, 

where \l is a positive number such that /i < min{l, n/h}. As before, we omit elementary calculations, 
noticing only that the leading terms in (|3.11|) can be obtained by the substitution in ( |3.7| ) instead of 



the function L the leading terms of its asymptotics for z — > +ioo. 



5. The solution a of equation ( |3.5|) is uniquely determined by its asymptotics (|3.10| ) - (|3.11|) , and 
by the fact that it is analytic and has no zeros in the strip Sq. 

Indeed, let o\ be another solution of ( |3~5| ) possessing these properties. By (|3.5|), the ratio ai(z)/a(z) 
a 2/i-periodic function. It is analytic in Sq, and thus, due to the periodicity, it is entire. On the other 
hand, it tends to 1 as z — > ±ioo. Therefore, this ratio is identically equal to 1. 

6. The cr-function satisfies the functional relations: 



17T 



(3.12) a(z + n) = (l + e h z ) a (z-n) 

(3.13) o-(-z) = e 4ft 2 + 12 + 12 



o (z 



a 



These three relations can be proved by using almost one and the same argument. To justify, for 
example, formula (|3.12|) , first, one checks that the ratio r (z) = a(z)/a(z — 2tt) is a 2/i-periodic 
entire function. Really, both the functions cr(z) and a(z — 2ir) satisfy (|3.5| ). This implies that r 
is 2/z-periodic. The function a(z) is analytic in the half-plane — 7r — h < Kez, and l/cr(z — 2tt) is 
analytic in the half-plane Re z < 3tt + h. So, r is analytic in the strip —it — h < Re z < 3tt + h. Thus, 
being 2/i-periodic, it is entire. Furthermore, the asymptotics of a imply that r (z) — > 1 as z — ► —too, 

and r (z) = e ~ mz l n + Z7r /"-(l + o (1)) as z — > -Moo. This is possible only if r (z) — 1 + e ft ft , 
which proves (|3.12|) . 



3.3. Meromorphic solution of equation (|3.4j) . Now, we come back to equation (|3.4j) . One of its 
solution can be constructed by the formula 

(3.14) v(p) = e L 2h P (T(p-po)/v(p + Po), 
where 

(3.15) Po = ii + h/2. 
Describe analytic properties of v. 



First of all, note that (|3.13| ) implies that 

(3.16) v (—p) = v (p). 

The function v has two chains of poles: 

±(-p + ir + h + 2nj + 2hk), j,k E NU{0}. 
If these two chains do not intersect, then the poles ±(— po + n + h) are simple, and 



( 3 17) »= J h c it 2 / 2h - ir£/ 2h - iir 2 / 12h + ih / 24 1 

1 ' 1 V 7T a{2i£ - ?r)' 

The asymptotic formulae for a imply 

(3.18) v(p) =e ±ipop/2h {l + o{l)), p^±ioo. 

3.4. Minimal solution of the model equation for Repo < tt- I n this section we assume that 

(3.19) Im£> -ir + h/2. 

This condition is equivalent to the inequality Rep < 71 ■ It implies that the poles of v are outside 
the strip 

(3.20) -h < Rep < h. 

1. Construct an entire solution of the auxiliary equation ( |3.2|) . Let Y be a vertical curve going 
along the imaginary axis from p = —iR.eC, to p = +iRe£ and having as the asymptotes for p — > ±ioo 

the line e~ m l A R. Since v has the asymptotics (|3.18|) , the integral in ( |3.3| ) converges and defines an 
entire function m (z). 



This function m is an entire solution of the auxiliary equation ( |3.2| ). To check this, we make the 
calculation from subsection |3.1|. As the convergence of the integrals is obvious, then, to justify this 



calculation, one has only to note that the horizontal distance from the contour Y to the poles of v is 
bigger than h. 

2. Let us turn the attention to the asymptotics of m for z — > ±ioo. One has 
Proposition 3.1. Let Im£ > — n + h/2. Then the function m has the asymptotics 
m (z) = a e^h ( z ~ n + + f (1 + (1)) + 

(3.21) +6 e~i( 2;_7r + ^) 2 + f (l + o(l)), z^+too, 

m(z) = c Q e2hi z - 7T - l S) -2(i + (l)) + 

(3.22) +e -2mz/h doe -l h (z-7T- i 2 -f (i + (i)), z ^ - ioo , 

where a , b , c and d are independent of z. The asymptotics are uniform in Rez if \Rez\ is 
bounded by a fixed constant. The asymptotic representations for w! can be obtained by differentiating 
the asymptotic representations for m. 

When proving the proposition, one checks also that the constant coefficients in the above asymptotics 
are given by the formulae: 

(3.23) a = 2iV^he~^h K ~ n ) 2 ~ K + !e, 

(3.24) b = 'J™ e-TH^-O 

a (2z£ - 7r) 

(3.25) c = -2v / ^e~^( 7T + ^) 2 ~^ + S, 



2Vvr/ii 2 < i U I if I inn 2 ih 



For the case where £ = 0, we have proved this proposition in |J. In the case under the consideration, 
one can use the same arguments and, even, almost the same estimates and calculations. So, we shall 
not repeat them here, and give only some comments on the proof. 

Begin with the case of z — ■> +ioo. Then, the asymptotics contains two leading terms. The first one 
is defined by the behavior of v as p — > — zoo, and the second one is related to the pole of the function 
v situated at po — it — h. Describe the first one. 

Remind that as p — » ±zoo, the integrand in (]3.3|) has the asymptotics 



_Hr» L r> 2 4- r> + n , , ^ 

e hP 4hP 2JiP(l+ (l)). 

If z — > +200, the integrand has a saddle point. It satisfies the equation 

i (-fp-4p 2 +ip-^p)=o, 

and thus, the saddle point is p\ = —2z + n — p . The first of the leading terms equals to the 



contribution of this saddle point to ( |3.3|) , e.i. to 

e~2k z% VAttHc'^ e~i zpi ~±P 2 i + ^~ ±P^ 



2 i iz 



= a eT h (z-* + ^ 2 + 2> 
where ao is given by ( |3.23|) . 

The second leading term of the asymptotics of m (z) equals to the contribution of the pole po — tt — h 
to (PQ, e.i. to 

r2 / A^lh „V,2 



e'TH z 2 2m Res p=Po ^ h (e~ iz Pl h ~ l P l Ah + m P/ 2h v (p) 



with b from (|3.24j ). 

The asymptotics ( |3.21|) of m (z) for z — > +zoo can be described as the sum of the contributions of 
the above two terms. 

In the case where z — > —ioo, the asymptotic ( |3.22| ) of m again contains two leading terms. The first 
one is due to the saddle point P2 defined by 

i {- l ip-±p 2 + f h p + ^p) = ^ 

i.e. P2 = —2z + n + p . The contribution of this saddle point to ( |3.3|) is given by 



e 



-k z2 V4^h e 3 -¥ e ~i m -IkP 2 2 + 2TPi + mPoP* 



where Co is given by ( j3.25|) . This is the first leading term, and the second one equals to 
-e~TH z 2 2m Res p= „ P0+7r+h (e~ iz Pl h ~ ip2 / Ah + ™ p / 2h v (p)) = 

= e - 2 ^/^ oe - 2 7r(----^) 2 -f. 

□ 

3. Recall that we can reconstruct the components of a vector solution of ( [3.1|) by the formule 
m-\{z) = m(z) and mo(z) = m(z — h). This vector solution is a minimal entire solution. Let us 



The theorems 1.1a and 1.1b imply existence of the canonical vector solutions /1.2 and 01,2 of ( |3.1| ). 
The asymptotics of their first components have the form 

(3.27) (/i, 2 )iW = e ± 2k^ + ^) 2 + 2 f + z-+*oc, 

0+ = *f - 7T, 

and 

(3.28) (^ 2 )iW = e ± 2k^ + ^) 2 - 2 f + z^-*oc, 

0_ = — if — 7T. 

Here, we have chosen some of the possible values of the parameters <p±. Comparing the asymptotics of 

) { 

vri2\z) 

is really a minimal solution of ( |3.1| ) corresponding to the chosen canonical bases. 
Below, we shall refer to m as to a minimal entire solution of d3.2|) . 

3.5. Minimal solution in the general case. In the previous section, we have constructed the 
minimal entire solution m(z) for f being in the half-plane ( |3.19| ). It is analytic in f in this half-plane. 
Now, we are going to continue it meromorphically on the whole complex plane of f and, in result, to 
prove 



Theorem 3.2. Equation ( |3.2|) has a solution m (2, f) entire in z and meromorphic in f . It has the 
following properties: 

(i) m is analytic in f in the upper half-plane of the complex plane f , and t/ie function a (2z£ — 
7r) m (2, f ) zs analytic in the half-plane Im£ < 0. 

fiij //a (2if — 7r) 7^ ; then the behavior of m for z — > ±?oo described by formulae ( |3.21|) - 
(|3.22| ) with coefficients a , b , c and d given by ( |3.23|) - (|3.26|) . 

Above, we have indicated explicitly the dependence of m on f . The rest of the section is devoted to 
the proof. 

1. Let us study analytic properties of m as a function of f . As we have already noted, it is analytic 
in f if Im£ > —71 + h/2. Continue meromorphically m(z, f) into the whole complex plane of f. 
The idea is to find a relation connecting the functions m (z, f), m (z, f + Z7r) and m (2, f + 2i7r) 
for Im£ > — 7r + h/2, and to use this relation for the continuation. Such a relation exists since all 
the three functions m (z, f), e' lnz ^ h m (z, f + ot) and m (2, £ + 2m) satisfy ( |3.2| ), and since the space 
of solutions of the equivalent to ( |3.2| ) equation ( |3.1| ) is a two dimensional module over the ring of 
/i-periodic functions. Prove 

Lemma 3.3. //Im£ > — n + h/2, then 

(3.29) a (2i£ - tt) m (2, = a (2if - 3tt) [a m (2, f + wr) + f3 m (z, £ + 2iic)], 

where a and f3 are entire functions of f and entire and h-periodic in z. 

Proof. As equation (|3.1|) and (|3.2|) are equivalent, then the space of solutions of (|3.2|) is also a two 



dimensional module over the ring of /i-periodic functions K, and, moreover, for any two solutions / 
and g of (|3.2|) , the expression 

{f(z),g(z)} = f(z + h)g(z) - f(z)g(z + h) 

is /i-periodic, and f and are linearly independent over K iff { f (z), q(z)\ 3^ 0. We call { f(z), q(z)\ 



Assuming additionally to the hypothesis of the lemma that £ satisfies the condition 

(3.30) l/o-(2z£-37r) ^ 0, 
e.i. that 

£ ^ -in + z/i/2 + «7rZ + ihk, I, k = 0, 1, 2, . . . , 

we check that e mz ' h m (z, £ + Z7r) and m (0, £ + 2in) form a base in the space of entire solutions of 
( |3.2j ). For this, we calculate their wronskian. We shall write / ~ g if / differs from g by an entire in 
£ factor having no zeros. 

First, we study the asymptotics of the wronskian for z — > ±zoo. They can be easily calculated by 
means of the asymptotics of m. Assuming that &o (£ + ^tt) an d ao(£ + 2«7r) are nonzero, we get 

(3.31) {e inz/h m{ Z) i + in), m (z, £ + 2m)} ~ 6 (£ + wr) a (£ + 2ivr) (1 + o (1)). 

2; — > +ZOO 

But, in fact, formula (|3.23|) shows that cto(0 7^ for all £ G C, and formula (|3.24j) implies that 
&o(£ + in) 7^ under the condition ( |3.30| ). This justifies (|3.31| ). Similarly, one shows that 

(3.32) {e iwz/h m(z, £ + in), m (2, £ + 2ztt)} ~ d (£ + itt) c (£ + 2in)(l + o (1)). 

z — > —ioo 

Since the wronskian is an entire /i-periodic function of z, (|3.31|) and (|3.32| ) imply that it is in fact 
independent of z, 

{e^ h m(z,£ + in), m (z, £ + 2m)} ~ d (£ + m) c (£ + 2m). 

Furthermore, taking into account ( |3.25| ) and ( |3.26| ), we see that 

(3.33) {e i7Tz/h m (z, £ + in), m(z,£ + 2in)} ~ 1/tr (2z£ - 3tt), 

and thus, under the condition (|3.30| ), e mz / h m(z, £ + in) and m (z, £ + 2i7r) are basis solutions. 
Therefore, in this case, m (z, £) is their linear combination, 

m (z, = « e i7r2/h m (2, £ + nr) + /3 m (2, £ + 2i?r) 

with the coefficients 

a = {m {z, f), £ + 2i7r)}/{e i7rz//l m(2, £ + z'tt), m{z, £ + 2m)}, 

= | e W^ m(z,£ + ztt), m (2, 0}/{e i7r * /ft m (2, £ + ztt), m (z, £ + 2m)}. 
One can calculate the wronskians in these formulae in the same manner as above, which gives 

2m 

{m (z, 0,m(z,£ + 2m)} = Wl (£) + w 2 (0 e h z , 

{e^ h m(z,^ + tn),m(z, S)} = w 3 (£), 

Wj ~l/a(2it-n), j = 1,2, 3. 

This and (|3.33| ) imply the desired relation (|3.29| ) under the condition ( |3.30|) . But, then it is valid for 
any £ in the half-plane Im£ > —n + h/2 as an equality of two meromorphic functions. 

□ 

2. Now, we note that the integral ( |3.3| ) remains analytic in £ and entire in z while the pole 
Po — 7T — h = i^ — 7T — h/2 of the function v remains in the left half-plane of the complex plane, and 
its pole — po + tt + h = — z£ + n + h/2 remains in the right half-plane. Thus, m is analytic in £ if 
Im£ > -tt- h/2. 

The above result implies that Lemma 3.3 remains valid for all £ in the half-plane Im£ > —n — h/2. 

3. Consider the function 



Since the poles of a (2z£ — tt) are in the half-plane Im£ > h/2, $ is analytic in £ in the strip 
— tt — /i/2 < Im£ < /i/2. Lemma 3.3 immediately implies 



Lemma 3.4. The function $ can be analytically continued in the half-plane Im£ < h/2. 

Proof. The right hand side in ( |3.29| ) is analytic in £ if —2tt — h/2 < Im£ < — tt + h/2. Therefore, $ 
is analytic in the strip — 2tt — h/2 < Im£ < h/2. By means of ( |3.12| ), one can rewrite formula ( |3.29p 
in the form 

$ (z, £) = a (z, £) $ £ + ztt) + /3 (z, £) (l + el & + 2m )) $ (z s £ + 2m), 

and this relation allows to continue $ analytically from the strip —2tt — h/2 < Im£ < h/2 in the 
whole half-plane Im£ < h/2. □ 

The lemma implies that 771(2;, £) can be meromorphically continued on the whole complex plane. 
We shall not distinguish between m and its meromorphic continuation. The m is analytic in £ in 
the upper half-plane of the complex plane, and the product a (2i£ — 7r)m(z, £) is analytic in the 
half-plane Im£ < 0. 

4. Remind that m is entire in 2; for Im£ > — 7r + h/2. As the meromorphic continuation of m into 
the whole plane of £ was obtained by ( |3.29| ), m is entire in z for all £ different from the zeros of 
cr(2i£ — tt), i.e. outside the poles of m. 

5. Since m(z + h, £),m(z, £) and m(z — h, £) are meromorphic in £, equation ( |3.2|) remains satisfied 
for all £ G C as a relation for meromorphic in £ functions. 

6. To complete the proof of the theorem, we have to check only its statement (ii). Consider 
the canonical Bloch solutions of ( |3.2|) described by ( |3.27|) - (|3.28|) . One can represent $ (z, £) = 
a (2z£ — tt) m (z, £) by their linear combination with ^-periodic coefficients. In a vicinity C + of +ioo, 

$ (z, = A (z, £) h(z, + B (z, £) / 2 (*, 0, 2 e C + . 

The coefficients A (.,£), -B (.,£) can be expressed in terms of the wronskians of $ and /i i2 : 

A = {$, / 2 }/{/i, /„}, B = {f\, / 2 }. 

Thus, they are analytic in £ in the same part of the complex plane as $. Also, all the Fourier 
coefficients of A and B are analytic in £ in the same domain as $, and, therefore are described 
by the same analytic formulae as for Im£ > — tt + h/2, so everywhere outside the poles of a. The 
asymptotics ( |3.21|) - ( |3.22|) show that 

A - a (£) a (2z£ - tt), 5 - 6 (£) a (2z£ - tt), 

Im£ > — 7r + /i/2, z — > +200, 

outside the poles of a (2i£ — 7r). So, the zero-th Fourier coefficients of A and B are equal respectively 
to ao<7 and &o0" for all £ outside the poles of <r(2z£ — 71"), and all the Fourier coefficients with the 
negative indices are zero on the whole complex plane of £. This implies the statement of the Lemma 
concerning the representation ( 3.21|) . The applicability of Q3.22| ) can be extended in the same way. □ 



3.6. Estimates of the solution m. Here, we just list some estimates for m (z) and m\z) immedi- 
ately following from the asymptotic representations (|3.21|) and (|3.22|) . These estimates will be used 
in the next section. One has 

(3.34) \m(z)\<CP(z), \m'(z)\ < C (1 + \y\) P (z), 

where 



where x = Kez, y = Imz and (p = Im£. These estimates are uniform in Rez if |Re^| is bounded by 
a fixed constant. 

4. Minimal solutions of a perturbed auxiliary equation 
In this section, we consider the equation 

(4.1) ip(z + h) +ip(z - h) + 2e^ cos zip (z) = w (z)ip (z), 

assuming that £ is complex number, and w is a meromorphic function satisfying the estimate 

(4.2) \w(z)\ <CeO-~f i )\y\, 0</i, \y\>Y, 
for sufficiently big Y. Here, as below, in this section, we use the notations 

x = Rez, y = Imz. 

Our aim is to prove the existence of a meromorphic solution of (|4.2| ) having the asymptotics repre- 
sentations of the same form as ones of the minimal solution m of equation (|3.2|) . 

Remind that m has poles if a(2i^ — tt) = 0, e.i. if 

2z£ = 2tt + h + 2hl + 27rm, I, m = 0, 1, 2, 

Let 

£ = {±(2tt + A + 2/iZ + 2vrm), Z, m = 0, 1, 2, . . . }. 
In this section, we assume that 

(4.3) Hi £ E. 



One can try to invert the operator being in the left-hand side ( |4.1|) . As we shall see, this leads to an 
integral equation of the form 

(4.4) V (z) = m (z) + I k (z, C)ip (C) d<Z, z e 7. 

J-y 

Here 7 is a vertical curve, and k is constructed in terms of m and m two linearly independent 
solutions of (|3.2| ). 



(4.5) n (», C) = (»■ ^ ™ ( f " ? f " W 1 «■ (0, 
where 

(4.6) ^ C) = ct gl^_i) +ij 

and {m, m} = m(^ + h)m(z) — m(z)fh(z + Zi) is the wronskian of m and rh. This integral equation is 
our main tool. Note that the kernel k can be considered as a difference analog of the resolvent kernel 
arising in the theory of differential equations: now, instead of the canonical ^-function, 8(y — 77), 
which equals to zero if y > rj and to 1 if y < r], one encounters the function 



First, we shall study the integral equation (|4.4j) , and then we shall study its solutions and, in 



particular, check that they satisfy Q4.1D . The analysis of the integral equation is quite similar to one 
we have carried out in the case of Harper equation ip (z + h) + ip (z — h) + 2 cos zip (z) = Eip (z), 
E = Const , in || . So, we describe in details only its new elements. 

The main results concerning (|4.1|) are formulated in subsection 4.4. 



4.1. Operator K. Now, we are going to describe precisely the construction of the integral operator 
from (14.41) (the solutions m and m, the integration contour), to describe the functional space and to 



4.1.1. The solutions m and in. Let to be the minimal solution of ( |3.2| ) constructed in Section 3, put 



(4.7) m{z, =m(z,0- 

Here, we indicated explicitly the dependence of to on £. Together with to, m is an entire in z solution 

of 

Remind that to is analytic in £ outside the set of zeros of o (2z£ — 7r). Thus, both to and to are 
analytic under the condition ( [4.31 ). 

The new solution has the asymptotics 

m (z) = (c (0 + o (1)) ^ ~ ^ + ^ + f + 



(4.8) +e!T ((i o (O + o(l))e2hl 2: - 7r + ^ +T 



and 



z — > +ioo, 



m (z) = (a (0 + o (1)) e _ 2ir ( z ~ n ~ *0 a ~ T+ 

if ~ ' f\ 2 iz 

(4.9) (5 (£) +o (l))e^^- 7r -^ "7, 

z — > — zoo. 

Here, do, &o, Co and do are the coefficients from the asymptotic formulae for to. 

4.1.2. The wronskian ofm and to. When proving Lemma 3.3, we have already discussed the space of 
solutions of equation fl3.2|), the linear independence and the wronskians of its solutions. We calculate 
the wronskian of to and to as we have calculated wronskians in the proof of Lemma 3.3. We recall 
that the wronskian is /i-periodic and get its asymptotics for z — > ±ioo (using the asymptotics of m 
and to). This leads to the result 

{to (z), m(z)} = e^coiO «o(0- 
This and formulae ( |3.23|) - ( |3.25| ) imply 

(4.10) {m(z), m(z)} = -Arrihe^ 2 . 

Note that this result means that the solutions m and m are linearly independent over the ring of 
/i-periodic in z functions. 

4.1.3. The integration contour. Let (p = Im£. The contour 7 is a strictly vertical curve which does 
not pass through any pole of w, comes from — ioo along the line 

(4.11) x = -<p, , 
and goes to +zoo along the line 

(4.12) X = 7T + ip. 



To motivate this specific choice of the asymptotes of 7, note that the estimates ( |3.34j ) show that, \m\ 
is minimal along 7: 

(4.13) \m(z)\ < Ce~\y\l 2 ^—, z 6 7, 

Po(z) 

(4-14) *>M = { e «\v\/h ,< . 

Note also that 

... _ L.I In . . 



4.1.4. The functional space. We define the operator K as the integral operator with the kernel Ql.5| 
acting in the space L2 (7, p), where p is the weight reflecting as the behavior of m and in along 7 
so the one of w, see 



(4.16) p(z) = e( l - ^\y\pl(z). 

4.1.5. Compactness of the operator K . The kernel of the operator K is quickly decaying along the 



contour 7: using the estimates ( |3.33| ) for m and m', one can show that 

(4.17) p 1/2 (z)\ K {z,C)\p- 1/2 (C)<C{l + \v\)e-2\y\-2\v\ i v = Im(, z, ( G 7. 

The estimate ( 4.17|) implies that the operator K is a Hilbert- Schmidt operator in £,2(7, p). 



4.2. Solutions of equation ((4.1|) in a vicinity of the curve 7. Note that, m G £2(7,2?). There- 



fore, at least, one of the equations 

(4.18) ip = m + Kip and ip = Kip, 

has a nontrivial solution ip G £2(7,;?). In the sequel, ip denotes this nontrivial solution. If both 
the equations have nontrivial solutions from Li2('y,p), then ip is one of them. Show that ip can be 
analytically continued into a vicinity of 7 up to a solution of equation QL1|) . 

4.2.1. Analytic continuation. One has 

Proposition 4.1. There exists a positive 5, such that ip can be analytically continued in the 
horizontal (h + 5) -vicinity 0/7. 

Proof. For ( G 7, the kernel k (z, () is analytic in z in the horizontal /i-vicinity Sh of 7. The 
integral, representing Kip, converges uniformly in z being in any compact subset of Sh- Therefore, 
the function Kip can be analytically continued in Sh- But then, ip, being a solution of ( [4.18[ ), can be 
also analytically continued in Sh- 

Let 5 be the horizontal distance from 7 to the closest to it pole of the function w. If 5 < h, then 
having proved that ip is analytic in the horizontal /i-vicinity of 7, one can deform the integration 
contour in the formula for Kip inside its horizontal 5- vicinity, and check that, in fact, ip can be 
analytically continued in the horizontal h + 5- vicinity. If 5 > h, then one deforms the contour inside 
its horizontal /i-vicinity to check that ip is analytic in 2/i-vicinity of 7 and so on. 

□ 

For brevity, we shall not distinguish ip and its analytic continuation. 

4.2.2. ip and equation ( |4.1| ). Note that the horizontal width of the horizontal (h + 5)-vicinity Sh+s 
of 7 is bigger than 2h. Show that ip satisfies equation (Op for z — h, z, z + h G Sh+s- 

Let (n f)(z) =f(z + h) + f(z-h) + 2et coszf(z). Then, H ip = 'HoKip, but 

{H K ip) (z) = 27riRes c=z K(z + h,()ip (£) = 
(m (z + h)rh (z) — m (z) rh (z + h)) 



{m, m} 

Thus, we come to 



w (z) ip (z) — w (z) ip (z). 



4.2.3. Asymptotics of ip for z — > +200. Let us describe the asymptotic behavior of ip in a vicinity of 
7. We begin by proving 

Proposition 4.3. Fix positive s% and e<i so that E\ < 2h and Ei < min{/i, 2ir/h}. In the horizontal 
Ei-vicinity 0/7, the solution ip admits the asymptotic representation 

(4.19) ip{z)=m{z){A + A 1 {z))+m{z){B + B 1 {z)), z -> +ioo, 
where 

(4.20) A = 8{ip) + — ^ J m(()w (£) ip (0 d(, B = --L J^m (() w (() ip (() d( , 

and the functions A\ and B\ satisfy the estimates 

(4.21) \A x (z)\, \B x {z)\ <Ce-^, y > Y, 
for sufficiently big Y , finally, W = {m, rh}, and 



(4.22) Sty) 



1, if ip — m-\- Kip, 
0, if ip = Kip. 



Remark. Note that the form of the above asymptotic representation depends on the integral 
equation for ip . 

The proof of this statement is almost the same as the proof of Proposition 3.6 from H. Here, we just 



outline the main idea of the proof. The function ip satisfies one of the equations from ( 4.18 ). Thus, 
on the contour 7, one can estimate it by estimating the right hand side of the equation. Using the 
inclusion ip G L 2 (7,p) and the estimate Q4.171) , one obtains 

\ip(z)\ < Const e'\ y/2 \ /p (z). 

The integrals 

r m(C)w (C) 1P (C) d(, [m(C)w (0 iP (C) 



quickly converge as £ —>■ +zoo, and one can obtain the asymptotics of Kip (z) for z —>■ +ioo inside 
the horizontal /i-vicinity of 7 by replacing in the kernel of K the function 9 (z, () with 2i. This leads 
to the answer 

if) (z) ~ Am(z) + Bfh (z). 
We omit the elementary estimates leading to Q4.21] ). To prove (|4.19|) outside the horizontal h- vicinity 



of 7, one has to write down a convenient representation for the analytic continuation of Kip outside 
7. For example, if z is to the right of 7, and the horizontal distance between z and 7 is between h 
and 2h, one can use the formula 



(K1P) (z) = 2m Res i=z _ h K (z, C) V> (0 + / « (z, C) V (0 ^C- 

□ 



7 



Using the representation ( |4.19| ) for ip and the asymptotics ( |3.21| ) -( |3.22| ) for m, one finally obtains 



the asymptotics of ip for z being in the horizontal ^-vicinity of 7: 
iP(z) = (a + o (1)) ek ( z ~ n + + f + 



L /'-v _ /rr _L it">2 1 1 _ 



with 

(4.24) a = Aa (O, b = Ab (£) + B^ffi, 

where do, &o and Co are the coefficients from the asymptotic representations (|3.21|) - (|3.22|) for m (z). 

4.2.4. Asymptotics of if) for z — > —ioo. Describing the behavior of if) for z — ioo, one proves first 
the representation 

2tti 

(4.25) = (C + Ci(^))m(z) + (D + D 1 (z))e ~ 2 m( 2 ), 26 5 £i (t), Imz^-oo. 

Here, C\ and are functions satisfying the estimates of the form (|4.21|) , and C and D are constant 
coefficients given by 

C = Sty), D = -JL J e +2 ^m (C) w (C) $ (C) d(. 

Remark. Let us emphasize that, now, the asymptotic representation strongly depends on the 
integral equation for if): if if) is a solution of the inhomogeneous equation, then C = 1; if if) is a 
solution of the homogeneous equation, then C = 0. 

The representation ( 4.25|) implies that, in the horizontal ei-vicinity of 7, 
ip (z) = (c + o (1)) ek ( z ~ * ~ l 2 - f + 

27T2 if ' t^\^ i 

(4.26) e^~ z (d + o(l))e~2h^ z ~ 7l ~^) ~2 Z , z -> -ioo, 

with 

c = Cc (O, d = Cd (O + Da (O. 
where a , c and do are the coefficients from the asymptotic representations Q3.21]) - (|3.22|) for m (z). 

4.3. Solution if) outside the vicinity of 7. 

4.3.1. Analyticity near ±ioo. Being analytic in the horizontal h + 5- vicinity of 7, the function if) can 
be meromorphically continued into the whole complex plane just by means of equation ( |4.1|) . Since 
w is analytic in some vicinities C± of ±ioo, if) is also analytic there. 

4.3.2. Asymptotics in C±. Discuss the asymptotics of if) as y — > ±00 for arbitrary fixed x. Begin 
with the case of z — > +ioc. As for equation ( |3.2j ), see subsection 3.4.3, one can construct for ( |4.1| ) 
the canonical basis solutions f\^ analytic in a vicinity C + of +ioo. This can be carried in the same 
way as in the proof of Theorem 1.1. One obtains 

(4.27) f 1>2 (z) = e ± i (* - * + + 5^(1 + (1)), y -> +ioc. 

In C + , one can represent if) by their linear combination with /i-periodic analytic coefficients 

(4.28) iP(z) = A(z)f 1 (z)+B(z)f 2 (z). 

Now, we recall that the asymptotic representation ( 4.23|) is valid in the ei-vicinity of 7 and that as e\ 
one can be choose any number from (0, 2h). In particular, this representation is valid in a horizontal 
vicinity of the horizontal length bigger then h, i.e. the period of A and B. Therefore, comparing this 
representation with ( [4.23j ), we see that the periodic coefficients A and B are bounded as z — > +ioo, 
and 

(4.29) A(z)^a, B (z) -> b, z -> +ioo. 

Now, as A and B are /i-periodic, (|4.27 ) - (|4.29| ) imply that the asymptotic representation fl4.23|) 
remains valid and uniform in Rez if Imz — >• +zoo and |Rez| is bounded by any fixed constant. 

Reasoning: in the same wav one proves that the asvmptotic representation (14.261) for ib as y — > 



4.4. Results of the section. The main results of this section can be formulated in the following 
form. 



Proposition 4.4. Let 2i£ ^ E. Let also w be a meromorphic function satisfying the estimate (Q|). 
Assume that 7 is a strictly vertical curve not passing through any pole ofw. Then equation ( }4.1| ) has 
a meromorphic solution ip {z) with the following properties: 

1. ip is analytic in z in some vicinities C± of±ioo; 

2. there is a positive 5 such that if) has no poles in the horizontal (h + S)-vicinity of r y; 

3. ip admits the asymptotic representations of the form 



i[> ( Z: f ) = ( a + o (1)) e&( z ~ n + l ^ + T + 



(4.30) 
and 

(4.31) 



+ (6 + o(l))e 2^ 



-4(z-7T + z0 2 + f^ 



y -> +00, 



i/;(z, = (c + o(l))e2ft 



4(z-7T-^) 2 -f 



2ni 

e h z (d + o(l))e 2ft 



2 + 



y 



—00, 



where a, b, c, and d are independent of z. The error estimates are uniform in x if \x\ is bounded by 
a constant. 



Remark. The solution ip satisfies one of the integral equations ( [4.18 ), if it satisfies the homogenous 
equation, then the coefficient c in the asymptotic representation ( f4.31| ) is zero. 

5. Minimal entire solutions of equation (|ITl"D 

In this section, we construct the minimal entire solutions of the matrix equation fll.lQ and prove 
Theorem 1.2. 

5.1. The plan. We start with the equation ( |1.5|) for the first component of a vector solution of 
(|1 . 1| ) . We represent in the form 



(5.1) 

where t satisfies the equation 
(5.2) 

This transforms (|1.5| ) to 

(5.3) 

with 

(5.4) 



^(z)=t(z)f (z), 
t(z + h) = p(z)t(z - h). 
f (z + h)+f (z-h)= Vl (z)f (z) 



t z) 



■v{z). 



t(z + h) 

We find a solution t of (|5.2|) such that t (z) / 1 (z + h) tends to some nonzero constants as z — > ±ioo. 
Then, we accept assumption ( |1.21| ). In this case, 



vi (z) ~ vi_ e mz + vi + e~ 



±2 OO, 



and, by a linear change of the variable, one transforms ( |5.3|) to equation (|4.1| ) investigated in the 
previous section. This allows to construct some analytic solution %pi of equation ( |1.5|) . This function 
can be considered as the first component of a vector solution ib of the original matrix equation (II. If). 



To be sure that ip is really entire in z, we check that it is analytic in a sufficiently wide horizontal 
vicinity a vertical curve 7. Really, since detM (z) = 1, the formulae 

(5.5) ip(z + h)=M(z)ip(z), if>(z-h) = M- 1 (z-h)i(>(z) 

show that, in this case, ip can be continued up to an entire function of z. 



5.2. Meromorphic solution of equation ( [5.21 ). Let 7 be a strictly vertical curve. Here, we 
construct a meromorphic solution of ( |5.2j ) analytic in a vicinity of 7. We shall use the following 
notations: 

• we denote by 7 + 2tt the curve obtained of 7 by the 27r-translation; 

• we denote by S y the strip bounded by 7 and 7 + 2n so that 7 C S" 7 and {7 + 2n} R 5 7 = 0; 

• we let n± = n±(6) and iV = n + + 

• we denote the zeros of b(z) situated in S y by z h I = 1,2, ... N. 

Proposition 5.1. Let 7 be a strictly vertical curve not passing through any point where either 



b(z) =0 or b(z — h) = 0. There exists a meromorphic solution t of equation (|5.2|) and a positive 
number 5 such that 

(i) t is analytic in the horizontal (h + 5) -vicinity 0/7; 

(ii) if Zq is a zero ofb(z) situated in this vicinity to the left 0/7, then t(z ) = 0; in addition, 
the multiplicity of the zero oft equals to the multiplicity of the zero ofb(z) at z = z ; 

(Hi) the function v\ defined by ( |5.4|) is analytic in the horizontal 8 -vicinity 0/7; 

(iv) the function t admits the asymptotic representations 

(5.6) t(z) =e in - z / 2 (l+T_(z)), Imz^-oo, 

(5.7) t(z) =t 00 e~ in + z l 2 {l + r + (z)), Imz^+oo, 
where 

(, n + +n_ 
% - J2 z i - ™N/2 + ihN/4 

and r± are functions satisfying the estimates 

(5.9) \t ± {z)\ < Ce-" |hz| , 

where \x is a positive number. This formulae are uniform in Kez if \Rez\ is bounded by a fixed 
constant. 



Proof. 1. Construct a solution of (|5.2|) in terms of the a-function introduced in subsection 3.2. 
The function p(z) = b (z)/b (z — h) can be represented in the form 

(5.10) p(z)=e m - h f[ Pl (z), Pl (z) l - e ~ <Z ~ Zl) 



I — e —i(z—zi—h) ' 
1=1 



Substituting this representation in equation ( |5.2j) , and comparing the result with ( |3.5|) , one immedi- 
ately finds out that the function 



is a solution of (|5.2j ). Clearly, it is meromorphic in z. 

2. To prove the proposition, we have to recall some properties of the cr-function. First, we note 



that its poles are at the points —n — h — 2hl — 27rm, I, m G N U {0}. Secondly, we remind that as 
the pole at z = —n — h so the zero at z = ir + h are simple. Moreover, if 2tti7i < 2h, m G N, then 
all the poles and zeros situated at the points z = —n — h — 2nl, and z = 7r + h + 2nl, I = 0, . . . m, 
are simple. 

The ratio has poles only at the points 

z = 2h + 2hl + 2-xm, and - 2-k - h - 2hl - 2-xm, l,m e NU{0}. 

The zeros of this ratio are at the points 

z = h + 2hl + 2Tcm, and - 2n - 2hl - 2-Km, l,m e NU{0}, 

the zeros at z — h and at z — —2-k are simple. Moreover, if 2irm < 2h, m G N, then all the zeros 
situated at the points z = h + 2irl and z = —2-k — 2nl, I = 0, . . . m, are simple. 

Now, we shall show that the described properties of the ratio imply the first three statements 

of the proposition. 

3. Show that there is a positive number 5 such that the constructed t is analytic in the horizontal 
h + 5- vicinity of 7. 

Denote by 8\ the horizontal distance between the set of the points Zi, I — 1, . . . , N, and the curve 
2n + 7. The poles of the factor ti(z) closest to 7 are at the points z\ + 2h and zi — 2ir — h. So, it is 
analytic in the horizontal h + 5-vicinity of 7, where 

5 = min{5i, h}. 

This implies the first statement of the proposition. 

4. Let us prove the second statement. Let z be a zero of b (z) situated to the left of the curve 7 
inside the above horizontal h + 5-vicinity of 7. Show that t (zq) = 0. 

All the factors t\ are analytic in the vicinity of 7. Inside this vicinity, to the left of 7, the factor t\ 
equals to zero only at the points z\ — 2nm, m = 1, 2, . . . . This proves the first part of the second 
statement. Prove the second part. 

The choice of 5 implies that a point of the form zi — 2irm, m = 2, 3, 4, ... , can be situated in the 
horizontal (h + 5)-vicinity of 7 only if 27rm < h. So, all the zeros of ti being in this vicinity to the 
left of 7 are simple. In result, the multiplicity of the zero of t (z) at any of these points being in the 
vicinity equals to the number of the factors t\ which are equal to there, e.i. to the multiplicity of 
the zero of b at this point. 

5. The function v\ (z) defined by ( |5.4j) is clearly meromorphic. It can be also represented in the form 

(5.12) „ lW = _^_ oW + _!<£)_ i(z - h) . 

Formula ( |5.12j ) shows that the poles of v± situated in the above horizontal vicinity of 7 can be only 
at the points where t (z + h) = or (and) t (z — h) — 0, e.i. at the points 

z = zi + 2hl + 2Tcm, and z = z x - 2-k + h - 2hl - 2irm, I,m6NU{0}. 

Among them, only the points 

z\ + 2nm and z\ — 2-nm + h, m = 0, 1, . . . , 

can be in the above vicinity. So, they can be only the points where either b (z) = or b (z — h) = 0. 
However, the curve 7 does not pass through any of these points. Thus, there is a positive 62 such that 
the function v-\ is analytic in the horizontal 5 9 -vicinitv of 7. If this 5? is smaller than the previously 



6. The asymptotics of t for z — > ±ioo follow immediately from formula (|5.11 ) and the asymptotics 



( |3.10| ) and ( |3.11| ) of the a-function. 

□ 

5.3. Entire solutions of ( |1 . 1| ) . Here, following the plan described above, we shall construct a set 
of entire solutions of 

1. Let 7 be a strictly vertical curve do not passing through any point where either b(z) = or 
b(z — h) = 0, and let t be the solution of equation ( |5.2| ) described in Proposition 5.1. Let us study 
in more details the function V\ defined by (|5.4j). 

By Proposition 5.1, v% is analytic in the <5-vicinity of 7. In the case of ( |1.21|) , the asymptotic 
representations for t imply that 

(5.13) Vl (z) = -e % ( nz + 0- + _ e -Knz + </>+ + tt) + w ^ 
where 

(5.14) </>_ = -n_(6) h/2 - zlnw_, 0+ = -n+(b) h/2 + i In v+, 
and w is a function satisfying the estimate 

\w(z)\ <Ce( n -ri\ lmz \, \lmz\ > Y, 



where \i is the same number as in ( |5.9|) , and Y is a sufficiently big positive number. This estimate 
is uniform in Re z if \Kez\ is bounded by a constant. 

2. Let 

(5.15) 0=(0 + + 0_)/2, 2^ = 0+-0- 
Equation ( |5.3|) can be rewritten in the form: 

(5.16) f Q (z + h) + f (z -h) + IS cos(nz + <\> + tt) f (z) = w (z) f (z). 
The change of the variable 

(5.17) z\ = nz + 4> + ir, 
transforms (|5.16 ) to the equation 



(5.18) i> {z\ + hi) + if} (zi - hi) + 2e ? cos(zi) ip (zt) = u>i(zi) i> (zi), 

where 

i> (zi) = fo(z), w 1 {z 1 )=w{z), h 1 =nh. 
This equation was investigated in the previous section. 

3. Now, we can construct a solution of ( |5.16| ) by means of Proposition 4.4. 



Check that the coefficients of (|5.18|) satisfy its assumptions. The estimate of w (z) implies that w\ 



satisfies an estimate of the form (|4.2|) as a function of Z\. Note that the branches of the logarithms 



in ( 5.14 ) are fixed modulo 2n. To ensure the assumption on £, we fix them so that the parameters 
4>± satisfy condition (|1.18|) . 



To apply the proposition, it rests to choose a vertical curve mentioned in it. We denote this curve 
71 and choose it so that it be the image under the transformation z 1— > z\ of the curve 7 used for 
constructing of the function t. To apply the proposition, we have to check that 71 does not not pass 
by any pole of w. But this follows from the statement (hi) of Proposition 5.1. 

Now, we construct the solution tb of (15.181) by means of Proposition 4.4. Note that this solution is 



determined by the horizontal distance between the poles of W\ and the curve 71. We can assume 
that 5i — nd ( if 5± > n5, we redefine it letting 5\ = n5; if 5i < n5, we redefine 5). 

Indicating explicitly the dependence of ip on £ and h, we let 

(5.19) f (z)=iji(noZ + <f> + n,Z,nh). 

The function /o is a solution of ( |5. 16 ). It is analytic in the horizontal h + 5- vicinity of 7. 



4. The function ^ = t f is a solution of equation ( |1.5| ), and, thus, can be considered as the 
first component of a vector solution \1/ of the matrix equation ( |1 . 1[ ) . Its second component can be 
reconstructed by formula Ql.6| ). 



5. Let us discuss the analytic properties of ipi- 

First, note that the asymptotics of fo (following from Proposition 4.3), and the asymptotics for t 
described by ( p.6| ) and (|5.7| ) imply the representations 

M*) = (Ao + o (1)) eit (nz + + ) 2 + f (n - n + ) + 

(5.20) +(B + o(l))e~2t( nz + ^ 2 + 'f ( n_n +), Imz^+oo, 
^(*) = (Co + o (1)) ("* + ^-) 2 " f (" - n -) + 

(5.21) +e -¥ z ( J D 1 + o(l))e"2^^ + ^-) 2 "f Imz^-oo, 
where n± = n±(b), and A , B , C , and Di are coefficients independent of z. 

Secondly, note that since both t and fo are analytic in the horizontal h + S- vicinity of 7, the function 
ip\ is also analytic there. 

6. Discuss the analyticity of the second component of the vector ^ = ^ ^ j . For any e > 0, 

we denote by L 7 (e) the part of the horizontal e- vicinity of 7 situated to the left of 7. Clearly, 
ipi(z + h) — 0(2)^1(2) is analytic in L 7 (/t + 5). So, if 6 has no zeros in this domain, ^2 (2) — 
(■01 (z + /i) — a (2;) i[)i(z))/b (z) is also analytic there. 

7. As the horizontal length of L^(h+S) is greater then h and $ is analytic here, then we can continue 
\1/ in the whole complex plane by means of equation ( |1 . 1|) : we continue \1/ to the right just by means 
this equation itself, and we continue it to the left by the formula ^(z — h) = M~ 1 (z — h)ty(z) (remind 
that detM = 1). 

We have proved 

Proposition 5.2. Assume that 6^0 and that v (z) satisfies condition (|1.21[) . By means of 
formulae (|5.14j) , choose some numbers <fr± satisfying ( |1.18|) . Let 7 be a strictly vertical curve not 
passing through any point where either b (z) = or b (z — h) = 0. If the curve can be drown so that 
there would be no zeros of the function b in the domain L 1 {h), then equation ( |1.1|) has an entire 
solution \l/ with the first component admitting the asymptotic representations (|5.20| ) - ( [5.21 ). 



In fact, the solution \I/ described in Proposition 5.2 is one of the main personages of this paper. In 
subsection 5.5, we shall see that it is a minimal entire solution of equation fll.lj) . 

5.4. Modified approach. In the previous subsection, we have constructed a (minimal) entire so- 
lution of equation( |l.l|) under a hypothesis on the geometry of the set of zeros of b(z). We were 
assuming that there is a strictly vertical curve 7 such that biz) and biz — h) have no zeros on 7 and 



this geometric condition. We begin with the case where b has only one simple zero zq G L 7 (h). Let 
us outline the idea. 



5.4.1. The idea. Any solution of ( |1.5|) can be considered as the first component ipi of a vector solution 
^ = ^ ^ j of equation (|1.1| ). Its second component can be recovered by ( |1.6|) . If ^1 is analytic in 
the horizontal (h + <5)-vicinity of a vertical curve 7 and such that 

(5.22) Vi (^0 + /i) - a (zo) ^1 (zb) = 0, 

then the corresponding vector solution \1/ is analytic at least in L^(h+5) and, so, it can be analytically 
continued up to an entire function just by means of 

In the case of Proposition 5.2, we have constructed (minimal) entire solutions of ( p. . 1|) in terms 
of a solution f of equation ( |5.3|) analytic in the horizontal (h + <5)-vicinity of 7, using the formula 
ipi(z) = t (z) fo(z). Now, by Proposition 5.1 (statement (ii)), 

(5.23) t(z ) = 0. 

So, = t fo appears to be analytic in the horizontal (h + <5)-vicinity of 7 even if fo is replaced by a 
solution of ( |5.3|) having there one (simple) pole situated at z . For the "old" analytic at z solution 
fo, the condition ( p\22| ) is not satisfied. The idea is to satisfy this condition by considering instead 



of the "old" fo linear combinations of the "old" fo and a "new" solution of ( |5.3|) having one simple 
pole at Zq- 

Of course, the meromorphic solutions of ( |5.3| ) we are looking for have to possed asymptotic rep- 
resentations for z — » ±200 of the same analytic structure as fo so that ipi would have the asymptotics 
of the form (|5T20| ) - QCTQ . 



To construct such solutions of ( p.3| ), we shall consider a solution m (z, zq) of the model equation 
1 possessing the analogous properties. Then, we shall plug it into the integral equation ( |4.4[ ) in 
the place of m (z) (the integral operator remains the same as before) and study its solutions. 

5.4.2. The solution m(z,zo). We shall use the variable z of the input matrix equation ( |1 . 1| ) . Let <fi 
and £ be as in ( [5.151 ). Consider the equation (compare it with ( pM6| )) 



(5.24) g(z + h) + g(z - h) + 2e ? cos(zi) g(z) = 0, z x = nz + (f) + 7r. 



The functions m(zi, £, hi) and rh(zi, ^, hi), hi = nh, are solutions of (|5.24|) . Here, we have indicated 
explicitly the dependence of these functions on h and £. 

Let 

m(zi) fh (zoi + hi) - m (z i + hi) m (zi) ] 



(5.25) m (z, z ) = 6 (z, z ) 



{m, fh} 



where {m, m} is the wronskian of m(z) and rh(z), see subsection |4.1.2j , 

z i = nzo + (j) + 7T, 

and, as before, 

6 (z, Zq) = ctg + 1. 

Since the wronskian {m, fh} is constant, and 6 (z, z ) is ^-periodic in z, the function m (z, z ) satisfies 



equation (|5.24|) 



Let us discuss analytic properties of m(z, zo). It is analytic outside the points z — z$ + Ih, I G Z. 
As [m(zoi) m ( z oi + ^1) — m (^oi + ^1) Th (201) ] = { m , ™}, the solution m(z, z ) has a simple pole at 

Z = Zo, 

h 



At the point z = zo + h, the numerator in (|5.25| ) equals to zero, and, thus, in fact, m(z, z ) is analytic 
here. So, for some positive 5, m (z, zq) is analytic in the horizontal h + 5-vicinity of 7 if z 7^ zq. 

The asymptotics of m(z) and fh(z) for z — > ±ioo imply the asymptotic representations for m(z, zq): 
m (z, z ) = A (z Q ) ( n2: + ^+) 2 + ¥ (1 + (1))+ 

(5.26) +5(z )e~^ (nz + ^+) 2 + 1 ir(l + o(l)), Imz^+ioo, 

m (z, z ) = e - 2 ™/ h (c (z ) ( nz + ^ ~ *f (1 + o (1))+ 

(5.27) + D(z )e~^ [nz + (l) -> ~~{l + o{\))\ , Imz^-zoo, 
where A, B, C and .D are independent of z. 

5.4.3. T/ie integral operator. To construct the meromorphic /o, we shall again use the integral op- 
erator from (|4.4| ), but now, it will be convenient to write it in terms of the variable of the input 
equation ( |1 . 1|) . So, now, 7 is a strictly vertical curve going to +ioo along the line Rez = — Re0 + / n, 
and from —ioo along the line Rez = — (tt + Re0_) / n,. We define the weight P(z) = e^ n ~^ y]{ Pq(zi), 



where y = Imz, /1 is the constant from the estimate for w (z) from ( |5.13| ), and po is defined by 
( 4.14Q . Both the functions m(.) and m (., Zq) belong to L 2 (7, P). 



The kernel of the integral operator takes the form 

/, „m / >n 1 fl/ >\ [ (gi) m (Ci) - m (Ci) m (zQ ] 

2m {m, mj 



Here, u> is the same as in ( |5.18|) . The kernel satisfies the estimate 



(5.29) P 1 / 2 (z)|/ ? (z,C)|P- 1/2 (C)<C(l + |r / |)e 2 If I 2 M y = Imz, V = lm(, 2, ( G 7. 
The integral operator .fT with the kernel ( |5.28| ) is compact in the space L 2 (7, P). 

5.4.4. Auxiliary construction. To construct the meromorfic /o, we use an operator equation slightly 
different from (|4.18| ). Before writing down this equation, we shall prove an auxiliary statement 
motivating this equation. Note that, since det M = 1, the equality b (z ) = implies that a (z ) 7^ 0. 
We put 

(5 30) s (z ) = - 1 — + ^ 

ft, f (z ) a (z ) ' 

One has 

Lemma 5.3. Assume that there exists a function f G L 2 ('~y, P) which can be analytically continued 
in a vicinity of the point Zq + h, and satisfies the relation 

(5.31) f(z) = Kf(z) + s(z )m(z, z )f(z + h)+m(z), z e 7, 

where either m (z) = m (zi) or m (z) = 0. Lei V ; i(- 2 ) = t i z ) f{ z )- This function has the following 
properties: 

(i) it is a solution of ( |1.5| ) analytic in the horizontal (h + 5) -vicinity 0/7, where 5 is a positive 
number; 

(ii) satisfies the relation (|5.22|); 



Remark. As we shall see from the proof, if the function f satisfying Q5.31J) exists, then it satisfies 
equation ( p.3| ) . It is the needed meromorphic solution of this equation. Furthermore, the function 
■01 = tf ^ the first component of the desired (minimal) entire vector solution of equation ( |1 . 1| ) . 



Proof, a. Since / G £2(7, P)i an d the kerenel of K satisfies the estimate ( |5.29| ), then the term Kf 
is analytic in the horizontal /i-vicinity of 7. Therefore, the right hand side of ( p.31| ) is analytic in 
this vicinity except the point z = zq. This implies that / itself is also analytic there. At the point 
z = Zq, this function has a simple pole, 

(5.32) res 2=20 f(z) = - s (z ) f(z + h). 



71 



Let Si be the horizontal distance from z to 7. Deforming the integration contour in the definition 
of K, one checks that, in fact, / is analytic in the horizontal (h + <5i)-vicinity of 7 without the point 

z = Zq. 

b. Using a calculation analogous to one from subsection [4.2.2] , by means of the residue theorem, we 
check that the relation ( |5.31| ) implies that, in the above vicinity of 7, / satisfies the equation (|5.16 ) 
which is just a different form of writing of (|5.3|). 



c. Remind that t is analytic in a horizontal (/i+5)-vicinity of 7, and t (zq) =0. If 8 > Si, we redefine 
it letting S — Si. Let ipi = t f. This function is analytic in the whole horizontal (h + <5)-vicinity of 7. 



Since / satisfies (|5.3|) , ipi is a solution of equation (jL 
d. Check that ipi satisfies the relation ( |5.22|) . Using (|5.32 ) and ( |5.30| ), we see that 

ipi(z ) = t' Oo) res z=20 f(z) = t _f° + A. J( Zq + h) = ipi{z + h). 

a (z ) a [zq) 



e. Finally, applying the methods of section 4 to equation ( 5.31 ), one obtains the asymptotics of / 
for z — > ±ioo which leads to (iii). 

□ 

5.4.5. Operator equation. Now, we turn to constructing a function / satisfying the hypothesis of 
Lemma 5.3. Since (|5.31| ) can not be used directly as an equation for a function from 1/2(7, P)i we 
change the functional space to Ti. = L 2 (7, P) ® C. For F G H, we use the notation 

/ 



, /GL 2 ( 7 , P), seC. 

Then we define the operator K, acting in Ti by the formula 

£ F= f K f m(z,z )s 

\ Kf(z +h) m (z + h, z ) s 

For / G L2 (7, P), the function Kf(z) is analytic in the horizontal h- vicinity of 7, and thus the 
number K f(zo + h) is well defined. 

One can easily check that the application / — ► K f(zo + h) is a linear bounded operator from 
L2 (7, P) to C, which implies that K, is a compact operator together with K. 

The new integral equation is 



(5.33) JC F + F = F, F (z) 



m [z) 
m (z + h) 



We let m(z) = m (zi) if 1 ^ spec fC and rh (z) = otherwise. 

Clearly, ( |5.33| ) always has a nontrivial solution. For F being a solution of this equation, 



(5.35) s = Kf (zq + h) + m (z + h, z ) s + m (zo + h). 

Discuss properties of the function /. 

5.4.6. The solution of the operator equation and the hypothesis of Lemma 5.3. Since / £ -^2(7, P), 
it can be analytically continued in a vicinity of zq + h just by means of ( 5.34Q . 

Letting in ( |5.34j ) z = z + h and comparing the result with ( |5.35| ), we see that 

s = f (zo + h). 

Substituting this expression for s in G5.34j ), we come to (|5.31|) . 



Thus, we see that / satisfies the assumptions of Lemma 5.3. This means that we have constructed 
the desired (minimal) entire solution. 

5.4.7. Several zeros ofb. We have proved the existence of the minimal solution in the case where b 
has only one simple zero in the domain L 1 (h). The case of several zeros (even multiple ones) can be 
treated similarly. In particular, if, in L^(h), there are J simple zeros of b, then instead of ( |5.31| ) one 
considers the relation 

j 

f(z) = Kf(z) + (zj) m (z, Zj) f(zj +h)+rh(z), z e 7, 
i=i 

and instead of Ti one has to introduce the space L 2 (7, P) © C J . We omit the details and formulate 
the result: 



Proposition 5.4. Assume that 6^0 and that v (z) satisfies condition Q1.21 ). By means of formulae 
( 5.14 ), choose some numbers <fi± satisfying ( |1.18| ). Equation (|1.1[) has an entire solution $ with the 
first component ipi admitting the asymptotic representations ( |5.20| ) - ( |5.21 ) . 



Remark. Note that iflE spec/C then, in the asymptotic representation ( |5.21| ) foripi, the coefficient 
C is zero 

5.5. Minimal entire solutions. 1. Now, let us discuss in more detail the behavior of the solution 
described in Proposition 5.4 for z — > ±ioo. We shall use the canonical basis solutions f\^ and 
analytic in some vicinities C + and C_ of +ioo and —zoo correspondingly, see Theorems la and lb. 



Remind that any entire solution ip of ( 1.1Q admits the representations (|1.19|) . Consider these repre- 
sentations for the solution \l/ described in Proposition 5.4. 

Remind that the asymptotics ( |1.10| ) of the basis Bloch solutions /i^ contain parameter <f> + . This 
parameter is defined by ( |1 . 1 1| ) modulo 2tt. The choice of the value of the parameter is a choice of 
two particular solutions fij. 

Formulae ( |1.11|) and (|5.14|) coincide, and thus, <p + from ( |1.10|) can be chosen equal to <p + in (|5.20| ). 
In this case, comparing ( |5.20|) and ( |1.10D , we see that in a vicinity of +ioo, admits the first of the 
representations ( |1.19| ) with A and B bounded as z — » +ioo. 

Choosing the canonical basis gi^ m a vicinity of — too with the same values of 0_ as in (|5.21 ), one 
also sees that ^ admits the second of the representations ( |1.19| ) with C and D bounded as z — > —ioo, 
and that, moreover, D (z) — > as z — >• —zoo. 

Now, we recall the definition of the minimal entire solutions of ( |1.1|) , and see that the solution \1/ is 
one of them. It is the solution ibjj. 



2. The proofs of the existence of the other three minimal solutions are similar to the above one. 



starting point for the analysis of Section 4. Consider for example the case of the solution for which 
B (z) — > as z —>■ +ioo. For this solution, equation (|4.4|) has to be replaced by 



(5.36) ip(z) = m(2ir-z) + / k (2tt - z, 2tt - C)V> (C) ^C- 

J{2tt-7} 

with the same kernel k as in Section 4. The contour {2ir — 7} is obtained of the contour 7 from 
Section 4 by the mapping z — > 27r — 2. As in Section 4, constructing a minimal solution of (|1 . 1|) , we 
have used as the first term in the right hand side of the integral equation the minimal solution of 
equation (|3.2| ) having the same form of the asymptotics for z — * ±ioo as the solution we construct. 
The analysis of equation Q5.36 ) is similar to one of equation (|47 



3. We have described how one can construct four minimal entire solutions corresponding to the 
canonical bases with the given parameters <j)±. The only condition on these to numbers is given by 
(|1.18p . In fact, it means that the canonical bases have to be consistent. This remark finishes the 
proof of our Theorem 1.2. 

6. Asymptotic coefficients of minimal entire solutions and the basic properties 
OF these solutions and of the corresponding monodromy matrices 

Choose two canonical bases f\^ and gi^- Consider the corresponding minimal solutions. Assume 
that these solutions exist and that their asymptotic coefficients are non zero. In this section, we study 
some basic properties of the minimal solutions and of the corresponding monodromy matrices. We 
shall use the notations from the section "Asymptotic coefficients" of the introduction. 

6.1. Wronskians of the minimal solutions. Here, we study the wronskians of the minimal entire 
solutions and check Proposition 1.3. One has 

Lemma 6.1. The wronskian of any two of the four minimal solutions corresponding to the bases 
/12 and gifl is independent of z. 

Proof. Consider the minimal solutions ip B and ip B . Study their wronskian. Remind that it is an 
/i-periodic entire function. 

In a vicinity of +ioo, 

{M*)> M*)) = {A {D \z)h{z) + B™{z) f 2 (z), A^ B \z) f x {z) + B^ B \z) f 2 (z)} = 

= w f (A^ (z) B^ B \z) - (z) B^ D \z)) . 

Here, we have marked by the letters B and D the functions A (z) and B (z) from the representations 
fll.l9f) for the solutions ipn and ips- The obtained formula for their wronskian implies that 

(6.1) {^d(z), iPb{z)} — > -w f A B B D , z —> +ioo. 
Similarly, one proves that 

(6.2) U>d(z), iPb(z)} — > w g C D D B , z — ► -zoo. 

Thus, the wronskian is a bounded entire function, and, therefore, it is independent of z. 

The cases of the other pairs of the minimal solutions can be treated in the same way. □ 



Asymptotics (|6.1| ) and ( |6.2|) imply the formulae 

(6.3) {ip D , ip B } = w g C D D B = -w f A B B D . 

In fact, this is two of the formulae described in Proposition 1.3. The other formulae can be derived 



6.2. Uniqueness of the minimal entire solutions. Here, we deduce from Proposition 1.3 its 
Corollary 1.5. It suffices to check 

Lemma 6.2. Consider two of the minimal solutions corresponding to the canonical bases fi i2 and 
gx t 2- If their wronskian is nonzero, then each of these solutions is unique up to an independent of z 
factor. 

Proof. Consider the case of the solutions ipo and if) B . The other cases can be treated similarly. 

Let if) be a minimal solution for which, as for ip£>, D (z) tends to zero as z — > —ioo. Show that 
if) (z) = Const ip D (z). By Lemma 6.1, the wronskian w = {iprj,ipB} is independent of z, and by the 
assumption, it is nonzero. Thus, if>£> and ip B form a base in the space of entire solutions of ( |1.1| ), 
and, so, if) can be represented in the form 

ip(z) = a (z) ip D + P {z)if)B, 
where the coefficients a and (3 are entire and /i-periodic, 

a(z) = - {ip (z),ip B (z)}, P(z) = - {if) D (z), V> (*)}■ 
w w 

Since all the three solutions are minimal solutions corresponding to one and the same pair of canonical 

bases, a and (3 are constant. Let C (z) and D (z) be the periodic coefficients from the representations 

(|i~20| ) for the solution ip, and let C (D) (z) and £> (D) (z) be the ones for the solution if>r>- As in the 

proof of Lemma 6.1, we check that 

(6.4) /*(*) = — (z) D{z)-C (z) D^ D \z)), 

w 

in a vicinity of —ioo. As D(z) and D^ D \z) tend to zero as z — >• —ioo, formula ( |6.4j ) implies that 
(3 — > as z — > —ioo, and, thus, (3 = 0. So, the minimal solution ipD is unique up to a constant factor. 
In the same way, one proves the same statement for ipg. 

□ 

6.3. Monodromy matrices corresponding to the minimal entire solutions. We begin this 
subsection by proving Theorem 1.6. 

Let 

w = {if) D , if> B }- 

By the hypothesis of the theorem the wronskian w is nonzero. It is given by formula d6[ 
The monodromy matrix corresponding to if>D and ips is defined by 

(V>2>(2 + 27t),VbCs + 27T)) = (lf) D (z),lf) B (z)) M T {z). 

The coefficients of the monodromy matrix admit the representations: 



(6.5) 



M u {z) = - {if) D (z + 2vr), if) B (z)}, M l2 {z) = - {if) D (z), ip D {z + 2tt)}, 
w w 



M 2 i{z) = - {ip B (z + 2vr), ip B (z)}, Mti{z) = - {if) D {z), ip B (z + 2tt)}. 
w w 



As in the proof of Lemma 6.1, one can easily calculate all the wronskians in ( |6.5|) using the canonical 
bases /i 2 and (?i >2 . Consider, for example, the coefficient M.\\. By means of the representation 
0, we get 

(6.6) M n = h ( A{D) (z + 2vr) A( B \z) {f\(z + 2vr), f\(z)}+ 

+A^ D \z + 2vr) B^ B \z) {h{z + 2vr), f 2 {z)} + B^ D \z + 2tt) A^ B \z) {f 2 {z + 2n), f 1 (z)}+ 



Now, we recall that fi >2 and gi >2 are Bloch solutions, 

f 1>2 (z + 2tt) = a lj2 (z) f lj2 (z), 9iA z + 2n ) = PiA z ) 9iA z )^ 
where ai )2 and j3\ y2 are /i-periodic functions. This allows to continue the calculation began in ( |6.6| ): 
M n = (ax(z) A&\z + 2vr) B^ B \z) - a 2 {z) B™{z + 2ir) A^ b \z)) . 

The last formula allows to get the asymptotics of M.\\ for z —> +ioo. Using the representation ( |1.16| ) 
and recalling that n+(v) = n, we get 

Mn = -<4 — A B B D e ft z (l + o(l)), z -> +ioo. 
w 

Here, A B and I?d are the asymptotic coefficients of ips and ip B . Now, using the formula for to, 
we get finally: 

(6.7) M n = a°e~^ 2 (l + o(l)), z -> +ioo. 



is nonzero. 



Remind that the constant a!] from ( |1.16 
Similarly, one can obtain 

1 27rin 

(6.8) Mn = ^e T_ T' (l + o(l)), — > —zoo, 

with the nonzero constant from ( |1.17| ). Formulae ( |6.7[ ) - ( |6.8| ) imply that the entire /i-periodic 
function A^n is a trigonometric polynomial. Considering it as function of the variable Z\ — ^r-z, we 
see that 

(6.9) n±(Mn) = n. 

Formulae ( |6.7D - ( |6.8| ) are, in fact, the first two of the formulae ( |1.25| ) - (|1.32j ) for the coefficients of 
the monodromy matrix corresponding to the solutions ipp and ipB- In the same way, one investigates 
all the other coefficients of the monodromy matrix. This leads to the formulae (|1.25|) - (|1.32| ). For 
M. considered as a function of z\ = z, they obviously imply that 

n_(M 12 ) < n - 1, n + (M 12 )<n, 

n-{M 2 i)<n, n+(M 2 i)<n-l, 

n±(M 22 ) < n — 1. 

These inequalities, and formulae (|6.9| ) mean that the monodromy matrix, as a function of the variable 
Zi — ^ z, belongs to Q (n). This completes the proof of the theorem. 

□ 

We finish this section by formulating one more theorem generalizing Theorem 1.6. As we have seen, 
any two of the minimal solutions form a basis in the space of entire solutions of ( |1 . 1| ) . We call these 
bases natural. One has 

Theorem 6.3. Let, the matrix M satisfy condition ( [1.21 ). For any natural basis, the coefficients 
of the corresponding monodromy matrix M. are trigonometric polynomials of the variable 

and for any of these polynomials the numbers n± satisfy the inequality 

—n <n±<n. 



The proof of this theorem is similar to the proof of Theorem 1.6. It can be said that Theorem 6.3 



6.4. Canonical factorizations of the monodromy matrices. For M G fi(n), n > 1, one can 
reconstruct in terms of the asymptotic coefficients of the minimal solutions only a part of the constant 
coefficients of the monodromy matrix. For example, formulae (|1.25 ) - (|1.26j ) allow to recover only 
the coefficients (Aiu) ±n of the trigonometric polynomial Ain, 

Mn(z) = (M 11 )_ n e ~ 2mnz / h + {M xl )_ n+1 e~ 2m ( n ~ 1 >/ h + 

+ ... + (M n ) n e 2mnz / h . 
Trying to get an efficient description of the other coefficients, we come to the following construction. 

Let /i : 2 and g± t 2 be two consistent canonical bases with the parameters <f>±. Consider canonical bases 
fil and gf\ with the parameters equal to <j>± + 27rj, j = 0, 1, 2, . . . n. Note that all the pairs f^\ 
and g^l are consistent. By definition, 

Jl,2 — Jl,2, 9l,2 — 9l,2, 

In view of the analysis of section 2.5.2, we can assume that 

/$(*) = 8+fo(* + 27 0> = s -9^( z + 27r )> 

where s± are two numbers, each of them can be equal either to +1 or to —1. For the sake of 
definitiness, we consider the case where s + — s_ — 1. 

Denote by ipp and ipg the minimal solutions corresponding to the canonical bases and g[^l- Let 

In view of Lemma 6.2, we can assume that 

^\z) =ij D (z + 2vr), ^\z) =^ B (z + 2tt). 

Let, for any j — 0, 1, 2, . . . n — 1, the solutions r/>^ and ^ be linearly independent over the ring 
of /i-periodic functions. 

Denote by the 2 x 2-matrix composed of the vectors and t/>£ . Any of the pairs ip% ] and ijjf 
is a basis in the space of entire solutions of ( |1 . 1|) . Therefore, one can define the transition matrices 
Sj(z) relating to 

*W(z) = *V- 1 \z)S t j (z), 

where "* " denotes the transposition. 

The monodromy matrix corresponding to the basis ipn, ipB is represented in the form 

M(z) = S n ^(z)S n . 2 (z) ... S^z). 

We call this factorization of the monodromy matrix canonical. The main feature of this factorization 
is related to the standard (canonical) form of the transition matrices Sj. One has 

Theorem 6.4. Any of the transitions matrices is an h-periodic trigonometric polynomial of the 
form 

aie 2niz/h + aQ + a _ ie ~^/h 0q + b _ ie -2niz/h 

c + C _ie 2 ^ 2 / h d 

where a±ip, 6_i,o, c^o and do are constant coefficients (depending on j), and a± ^ 0. The det Sj are 
nonzero and independent of z. 

The theorem shows that anv of the transition matrices coincides up to a constant factor with a matrix 



Sj can be expressed in terms of the asymptotic coefficients of the minimal solutions if) 



i) 



i) 



(?) 

D i 



4$, 



and 



by formulae similar to (|1.25|) - (P-.32). 



7. One dimensional difference Schrodinger equations 
Consider the difference Schrodinger equation 

where v is a given trigonometric polynomial, and E is a spectral parameter. It is equivalent to 
equation (\1 . 1 ) with the matrix 

(7.2) M=( 2£ - 1 2 "W - 1 

Assume that n + (v) — n_(f ) = n G N. Then, this matrix belongs to the set Vt nn . Thus, one can apply 
Theorem 1.6 to investigate the spectrum of ( |7.1[) by means of the monodromization procedure. 

In this section, we concentrate on the Harper equation 
(rj q\ f (x + h) + f (x - h) 

(7.3) + Xcoszf(z) = Ef(z), 

in which, A is a fixed positive parameter. Let us study for this equation the monodromy matrix 
described in Theorem 1.6 in more detail. 



7.1. Monodromy matrices for Harper equation. 1. First, discuss the choice of the bases fx )2 
and gi t 2- The formulae ( |1 . 1 1[ ) and ( |1.14j ) allow to choose the corresponding parameters <f>± so that 

(7.4) 0_ = -i£-7r, + = ^-vr, £ = In A G K. 

Since (el, these canonical bases are consistent, and there exist all the four minimal entire solutions. 

2. Consider the minimal solutions ipD and ipB corresponding to these canonical bases. Assume that 
these solutions are linearly independent. 

The solutions ip£> and ipB are defined up to independent of z factors. Describe a convenient choice 
of these factors. Since the wronskian {ipD, iPb} 7^ 0, the asymptotic coefficient Co is nonzero, see 
Proposition 1.3. Normalize the solution ip£> by dividing it by its asymptotic coefficient Cd- We use 
for the new minimal solution and its asymptotic coefficients the old notations. In particular, now, 
C D = 1. 

3. To fix uniquely the second minimal solution, we note that the matrix ( [7.2| ) with v = cos z satisfies 
the relation 

(7.5) M(2n-z) = aM- l (z)a, a = 

This implies that the vector 

(7.6) if>(z) = ai) D {2-n - z + h) 

satisfies ( |1.1|) together with ipr>- Show that, up to a constant factor, if) coincides with if>B- Begin with 

Lemma 7.1 For z being in a vicinity of — ioo, 

af^n-z + h) = -g 1 {z)c x {z), a / 2 (2tt - z + h) = -g 2 {z)e 2 {z), 




Proof. Let us prove the first of the above relations. Note that, for the matrix (|7.2|) corresponding 
to Harper equation, the general asymptotic formulae for f\ and g\ take the form 

{(■() hW-e ^ -A e"** (l + o(l)) 

Im z — ► +00, 

and 



(7.8) <7x(z) = e" 



_l e -i^(l + (l)) 



Im 2 — > — 00. 



Substituting ( |7.7| ) in the expression —a /i(27r — 2 + /1), we check that it does have the representation 
(|7.8j). Moreover, one can easily see that if f\ is a Bloch solution, then cr fi(2n — z + h) is also a Bloch 



solution of fll.ll) . Now, Lemma 2.6 implies the desired result. The second relation can be derived 
similarly. 

□ 

This lemma immediately leads to 

Lemma 7.2. The solution coincide with ip B up to a constant factor. 

Proof. The representations Q1.1SQ for ip B imply that, in a vicinity C_ of —ioo, 

(z) = A (d) (2tt -z + h)a / x (27r - z + h) + B {D) (2n - z + h) a f 2 (2ir -z + h). 
Therefore, in view if Lemma 7.1, 

(7.9) 4>(z)=C(z)g 1 (z)+D{z)g 2 (z), 
where 

C(z) = -ci(z) A (D) (2ir - z + h), D (z) = -c 2 (z) B 1 - d \2tx - z + h). 
Similarly, one shows that in a vicinity C + of +ioo, 

iP(z) = A(z)f 1 (z) + B(z)f 2 (z), 

where 

(7.10) A (z) = —a[ \2u - z + h)C {D) (2Tr - z + h), B (z) = — c 2 1 (2n- z + h) D {D) (27r -z + h). 

These representations show that tp is a minimal solution for which B (z) —>■ as z — > +zoo. Therefore, 
by Corollary 1.5, it coincides with ips up to a constant factor. 

□ 

4. In the sequel, we shall normalize the solution ipg so that 

(7.11) Tj) B (z) = aiJ)D(2Tr-z + h). 
Note that, in this case, ( |7.9|) and (|7.10j) imply the relations 



(7.12) A B = -C D = -1, B B = -D D e" 4 - 2 ^, C B = -A D , D B = -B D . 

5. Discuss now the monodromy matrix corresponding the chosen if) B and ip B . 

Theorem 7.3. The constructed minimal entire solutions ipo andip B and their asymptotic coefficients 
are meromorphic in E, 



(7.13) 

where 



M 



±st 

Ai 



and the corresponding monodromy matrix has the form 

( a-2X lC os(27rz/h) s + t e - 2mz / h \ 

K - s ^ te 2mz/h 

x l-s 2 -t 2 
a = Al- 
and 



st 



X 1 = X 2n/h , 



s = -Xi— -, t = -XiA D . 

Proof. First, prove that ipr, and ips and their asymptotic coefficients are meromorphic in E. Consider 
the solution ipD before its normalization by the condition Cp = 1. It can be constructed by means 
of the integral equation (|4.4p . For the matrix (|7.2|) with v (z) = cosz, this equation takes the form 



ij)(z) =m (z) + EKtp (z) 
where K is the integral operator with the kernel 

in {m, m\ 

This equation describes the first component of the vector The integral operator is compact in 
the suitable function space, and the solution iffp together with ijj appears to be meromorphic in E. 

Furthemore, the solutions f± t 2 and gi^, being constructed as in Section 2, are entire in E. 

Finally, the asymptotic coefficients of ip can be expressed in terms of the wronskians of the canonical 
basis solutions fi^ and and of ip itself, for example, C is the zeroth Fourier coefficient of the 

function 



{^ 1) ,92} 

{91=92} 



. Thus, the asymptotic coefficients are meromorphic in E. 



Note that we have already calculated these asymptotic coefficients for E = in Section 3: formulae 
( 3.23Q - ( p. 26 ) imply that these meromorphic functions are not identically zero. 

All this shows that the normalized minimal solution i/jp, the minimal solution ips related to it by 
(|7.11 ), and the asymptotic coefficients of these solutions are meromorphic in E. Moreover, since the 
asymptotic coefficients oiipo are nontrivial meromorphic functions, the wronskian {ipD, ?Pb} is not 
identically zero. 



Now obtain the representation (|7. 13|) for the monodromy matrix corresponding to the chosen minimal 
solutions ipD and ipB- First, consider the coefficient M.\\. By Theorem 1.6, it is a trigonometric 
polynomial of the form 

a _ ie -^iz/h + ao + aie 2mz/h 



Formulae flOo] ) - QEZg ) imply that 



a_ 1= a 2 , a 1 = p 1 , 

where a° and f}® are the first Fourier coefficients of the Bloch multipliers a<i and f3\ of the Bloch 
solutions f 2 and g 1: 

f 2 (z + 2tt) = a 2 {z) Hz), 9l {z + 2vr) = ^{z) 9l {z). 

The solution g\ admits the representation ( |7.8|) which implies that fi\ = —X 2n /^ = —X\. Similarly, 
a 2 = — Ai, and, thus, 



Formulae ( |1.27| ) - Ql.32j ) and relations (|7.12|) imply the desired representations for .M12, and 
M.22- Formula for the coefficient a® follows from these representations and the equality det.M 
I . □ 



7.2. Family of matrices generated by Harper equation. Fix a positive A. Define a family 
H(A) C Q (1) of matrix functions. This family consists of a two dimensional manifold H° and four 
one dimensional linear manifolds and hf. The H° is the set of the matrices -M(., A, w) defined 
by 

(a - 2A cos z s + te~ lz 
-s -te iz i st 

where 

1 - s 2 - t 2 
a = A - 



st 

w = (s,t) G C 2 , s,t^0. 



The Iiq is the set of the matrices 




M±(z, A, a) = I I , a e C, 



and hf is the set of the matrices 

(a — 2A cos z ±e~ 1 ' 
I , a E C. 
Te iz 

Apply Theorem 1.6 to equation ( |1.1|) with a matrix Mel Choose the bases /i^ and gi^ so that 

<f>- = -i£-n-~n-(b), </>+ = i£ - ir - - n + {b), f = lnAeK. 

These canonical bases are consistent, and there exist all the four minimal entire solutions. Assume 
that the asymptotic coefficient Co is nonzero. 

We normalize the solution ipjj by the condition (7/5 = 1 and define ipB by formula ( |7.11|) . Since the 
matrix M satisfies the relation (|7.5|) , ipB is really the desired minimal solution. 

Assume that ipo and ipB are linearly independent over the ring of /i-periodic functions. Instead of 
Theorem 7.3, one can prove 

Theorem 7.4. As a function of the variable 2-rrz/h, the monodromy matrix M. corresponding to the 
constructed minimal entire solutions belongs to H(Ai) with 

X, = \ 2 ^ h . 
In particular, if Ai 6 H° , then it has the form 

M{2irz/h, Ai, Wi), w 1 = (s 1 ,t 1 ), 
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